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Abstract  The  accuracy  of  state-of-the-art  global  barotropic  tide  models  is  assessed  using  bottom 
pressure  data,  coastal  tide  gauges,  satellite  altimetry,  various  geodetic  data  on  Antarctic  ice  shelves,  and 
independent  tracked  satellite  orbit  perturbations.  Tide  models  under  review  include  empirical,  purely 
hydrodynamic  ("forward"),  and  assimilative  dynamical,  i.e.,  constrained  by  observations.  Ten  dominant  tidal 
constituents  in  the  diurnal,  semidiurnal,  and  quarter-diurnal  bands  are  considered.  Since  the  last  major 
model  comparison  project  in  1997,  models  have  improved  markedly,  especially  in  shallow-water  regions 
and  also  in  the  deep  ocean.  The  root-sum-square  differences  between  tide  observations  and  the  best 
models  for  eight  major  constituents  are  approximately  0.9,  5.0,  and  6.5  cm  for  pelagic,  shelf,  and  coastal 
conditions,  respectively.  Large  intermodel  discrepancies  occur  in  high  latitudes,  but  testing  in  those  regions 
is  impeded  by  the  paucity  of  high-quality  in  situ  tide  records.  Long-wavelength  components  of  models 
tested  by  analyzing  satellite  laser  ranging  measurements  suggest  that  several  models  are  comparably 
accurate  for  use  in  precise  orbit  determination,  but  analyses  of  GRACE  intersatellite  ranging  data  show  that 
all  models  are  still  imperfect  on  basin  and  subbasin  scales,  especially  near  Antarctica.  For  the  M2  constituent, 
errors  in  purely  hydrodynamic  models  are  now  almost  comparable  to  the  1980-era  Schwiderski  empirical 
solution,  indicating  marked  advancement  in  dynamical  modeling.  Assessing  model  accuracy  using  tidal 
currents  remains  problematic  owing  to  uncertainties  in  in  situ  current  meter  estimates  and  the  inability  to 
isolate  the  barotropic  mode.  Velocity  tests  against  both  acoustic  tomography  and  current  meters  do  confirm 
that  assimilative  models  perform  better  than  purely  hydrodynamic  models. 
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1.  Introduction 

Tidal  signals  permeate  the  entire  Earth  system.  An  accurate  knowledge  of  tides  is  consequently  needed  for 
an  astonishing  variety  of  geophysical  fields  ranging  from  the  orbit  of  the  moon  to  the  mixing  of  the  oceans, 
from  solid-earth  geophysics  to  coastal  flooding  [e.g.,  Lambeck,  1 988;  Munk,  1 997;  Wilhelm  et  al.,  1 997; 
Vlasenko  et  al.,  2005],  An  especially  important  application  for  accurate  tide  models  is  providing  tide  "correc¬ 
tions"  to  various  measurements  so  that  smaller  nontidal  signals  may  be  studied.  For  example,  barotropic 
tide  models  are  used  regularly  to  remove  tidal  variability  from  space  geodetic  observations;  this  is  a  critical 
necessity  for  successful  satellite  altimetry  [e.g.,  Fu  and  Cazenave,  2001]  and  satellite  gravimetry  [Seeber,  2003; 
Visseretal.,  2010],  and  in  both  cases  improved  tidal  corrections  lead  to  a  reduction  of  aliased  tidal  "noise"  in 
nontidal  signals  of  interest. 
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Much  progress  has  been  achieved  in  recent  years  in  improving  global  tide  models.  In  no  small  part,  this 
progress  has  been  both  motivated  by  and  the  result  of  satellite  altimetry.  For  more  than  two  decades  altime¬ 
try  has  been  providing  critical  information  needed  to  constrain  models.  Yet  substantial  progress  is  still 
required  to  reduce  errors  in  presently  used  models.  For  example,  errors  in  basin-scale  tide  prediction  are 
roughly  comparable  to  measurement  errors  of  present-day  satellite  gravimetry  and  they  exceed  expected 
errors  of  future  gravimeters  by  an  order  of  magnitude  or  more  [ Koop  and  Rummel,  2008;  Wiese,  201 1;  Wiese 
etal.,  201 1],  Other  model  inadequacies  will  become  evident  in  the  discussions  below. 

Since  the  pioneering  works  of  Pekeris  and  Accad  [1 969],  Schwiderski  [1 979],  and  Parke  and  Hendershott 
[1980],  several  generations  of  tide  models  have  passed,  each  leading  to  further  improvements.  The  last  truly 
comprehensive  assessment  of  global  ocean  tide  models  was  published  by  Shum  et  at.  [1 997,  henceforth 
Shum97]  only  a  few  years  after  the  launch  of  the  Topex/Poseidon  (T/P)  satellite,  a  mission  whose  orbit  was 
specifically  designed  to  measure  tides  —  which  it  did  with  remarkable  success.  Shum97  documented  the 
degree  of  improvement  of  global  tide  models  that  resulted  from  the  analysis  and  assimilation  of  only  a  few 
years  of  T/P  data,  aiming  in  particular  to  identify  the  next  model  to  be  included  in  the  altimetric  Geophysical 
Data  Record  (GDR)  processing.  For  that  purpose  they  used  a  variety  of  tests  based  on  comparisons  to  in  situ 
data,  including  tide  gauges  and  bottom  pressure  recorders,  as  well  as  superconducting  gravimeters  and 
satellite  measurements. 

During  the  1 7 years  since  Shum97  was  published,  the  field  has  advanced  considerably.  We  have  witnessed 
not  only  improved  global  tidal  atlases  but  also  a  much  deeper  insight  into  the  dynamics  of  global  tides. 

Our  understanding  of  energy  dissipation,  the  role  of  baroclinic  processes  in  the  energy  budget,  and  the 
subsequent  fate  of  tidal  energy  and  its  wider  implications  have  all  seen  marked  progress,  and  these  topics 
continue  to  be  at  the  forefront  of  modern  research. 

The  central  goal  of  this  paper  is  to  update  the  review  of  Shum97  and  document  the  considerable  advances 
in  global  barotropic  tidal  models  themselves.  All  global  models  analyzed  here  were  submitted  by  develop¬ 
ers  willing  to  have  their  work  subjected  to  extensive  tests  and  comparisons.  The  main  focus  of  the  paper 
lies  in  the  sea  level  representation  of  the  main  tidal  constituents  in  modern  tide  models,  including  M2,  S2, 
N2,  K2,  01,P„Q1,K1,  M4,  and  MS4  (see  standard  texts,  e.g.,  Pugh  [1987]  and  Pugh  and  Woodworth  [2014], 
for  the  periods,  relative  amplitudes,  and  astronomical  origins  of  these  constituents).  The  primary  models 
being  tested  are  all  barotropic  models  tightly  constrained  by  satellite  altimeter  data,  and  these  models  are 
all  found  to  be  very  accurate,  even  by  standards  of  only  a  few  years  ago.  We  employ  many  of  the  same  tests 
performed  by  Shum97,  but  extensive  improvements  have  been  made  to  some  of  the  test  data.  For  example, 
the  expanded  tsunami  warning  network  [Bernard  etal.,  2010]  has  provided  new  seafloor  pressure  measure¬ 
ments  in  regions  previously  lacking  in  situ  data.  New  data  sets  are  now  available  over  polar  seas,  including 
valuable  new  tidal  GPS  measurements  obtained  on  floating  ice  shelves  [King  etal.,  201  la,  201 1  b]  as  well  as 
new  satellites  capable  of  providing  critical  polar  observations,  such  as  the  ICESat  and  CryoSat-2  altimeter 
missions  and  the  Gravity  Recovery  and  Climate  Experiment  (GRACE)  gravimeter  mission. 

With  these  new  test  data  sets  in  hand,  we  took  the  opportunity  to  test  several  recent  purely  hydrodynamic 
numerical  tide  models  (hereafter  also  referred  to  as  "forward,"  "unconstrained,"  and  "non-data-assimilating" 
models),  which  are  unconstrained  by  any  tidal  observations.  Comparisons  of  the  unconstrained  hydro- 
dynamic  models  are  limited  to  the  M2  component.  Finally,  we  attempt  an  assessment  of  the  accuracy  of 
barotropic  tidal  currents,  although  not  all  models  include  estimates  for  currents.  This  assessment  builds  on 
earlier  efforts  by  Luyten  and  Stommel  [1991],  Dushaw  et  at.  [1997],  and  Ray  [2001]  but  is  limited  mostly  to 
the  M2  constituent.  The  sections  on  forward  models  and  tidal  currents  include  results  from  two  recent  mod¬ 
els  which  simultaneously  simulate  tides  and  the  oceanic  general  circulation;  although  such  models  include 
internal  (baroclinic)  tides,  the  evaluations  of  such  models  here  are  confined  to  barotropic  elevations  and 
currents,  consistent  with  the  emphasis  of  this  paper. 

One  important  subject  not  addressed  here  is  the  long-period  tides,  which  are  forced  by  the  zonal  potential 
with  periods  from  1  week  to  1 8.6  years  [Wunsch,  1 967].  In  general,  these  constituents  are  smaller  than  the 
short-period  constituents  considered  in  this  paper,  and  fewer  models  have  been  developed  and  released. 
The  main  reason  for  the  omission  of  long-period  tides,  however,  is  the  lack  of  high-quality  test  data  sets  of 
harmonic  constants  comparable  to  those  available  for  short-period  tides.  In  particular,  in  the  deep  ocean, 
the  most  reliable  test  data  set  is  probably  still  that  developed  by  D.  Luther,  based  on  the  analysis  of  24  Pacific 
island  tide  gauges  [Miller  et  at.,  1 993].  There  is  great  need  for  a  more  comprehensive  global  test  data  set  for 
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Table  1.  Participating  Ocean  Tide  Models 

Model  Type3  Resolution  Authors 

Modern  data-constrained  models 
GOT4.8  E  1/2° 

OSU12  E  1/4° 

DTU10  E  1/8° 

EOTIIa  E  1/8° 

HAM12  H  1/8° 

FES12  H  1/16° 

TPX08  H  1/30° 


Historical  models 

NSWC 

H 

1° 

Schwiderski  [1 979] 

CSR3.0 

E 

1° 

Eanes  and  Bettadpur  [1 996] 

Purely  hydrodynamic  models 

HIM 

N 

1/8° 

Arbic  et  al.  [2008] 

OTIS-GN 

N 

1/8° 

Green  and  Nycander  [201 3] 

STORMTIDE 

N 

1/10° 

Muller  et  al.  [2012] 

OTIS-ERB 

N 

1/12° 

Egbert  et  al.  [2004] 

STM- IB 

N 

1/12° 

Hill  etal.  [2011] 

HYCOM 

N 

1/12.5° 

Arbic  et  al.  [201 0,  updated] 

aE,  empirical  adjustment  to  an  adopted  prior  model;  H,  assimilation  into  a  barotropic  hydro- 
dynamic  model;  N,  purely  hydrodynamic  (no  data  constraints). 


Ray  [1 999,  updated] 

Fok  [2012] 

Cheng  and  Andersen  [2011] 
Savcenko  and  Bosch  [2012] 
Taguchi  et  al.  [2014] 

Lyard  et  al.  [2006,  updated] 
Egbert  and  Erofeeva  [2002,  updated] 


long-period  tides,  but  this  requires  very  long  time  series  to  extract  these  small  signals  from  the  red-noise 
background  variability. 

The  remainder  of  the  paper  is  structured  as  follows:  in  section  2  we  briefly  summarize  all  models  involved  in 
the  quality  assessment.  Section  3  provides  a  summary  assessment  of  differences  between  data-constrained 
(empirical  and  hydrodynamic)  tide  models.  Section  4  evaluates  the  models  against  in  situ  data,  separately 
for  the  nonpolar  deep  ocean  and  for  shelf  regions.  Except  for  a  few  stations  involving  one  model,  these  in 
situ  data  are  independent  of  all  tested  models.  Section  5  focuses  on  an  evaluation  of  tide  models  in  high 
latitudes  based  on  tide  gauge,  bottom  pressure,  and  GPS  data.  Section  6  is  concerned  with  testing  the  tide 
models  against  various  kinds  of  satellite  data,  including  satellite  gravity  data,  which  are  completely  indepen¬ 
dent  of  all  models,  and  satellite  altimeter  data,  some  independent  of  the  models  and  some  not.  Section  7  is 
concerned  with  purely  hydrodynamic  models  and  their  ability  to  simulate  observed  tidal  heights.  Section  8 
then  tests  barotropic  tidal  currents  from  data-constrained  and  forward  models  against  estimates  avail¬ 
able  from  moored  current  meters  and  from  acoustic  tomography.  A  summary  and  an  outlook  are  provided 
in  section  9. 

2.  Summary  of  Participating  Tide  Models 

The  following  subsections  give  a  brief  overview  of  the  tide  models  evaluated  in  this  study.  We  structure  the 
discussion  into  (i)  seven  modern  (hydrodynamic  and  empirical)  tide  models  involving  data  (assimilation), 

(ii)  two  historical  tide  models,  and  (iii)  six  hydrodynamic  forward  models  simulating  tides  by  imposing  the 
respective  forcing  without  any  data  constraints. 

2.1.  Modern  Tide  Models 

Modern  tidal  models  can  be  classified  as  (1)  empirical  adjustment  to  an  adopted  prior  model  and  (2) 
barotropic  hydrodynamic  models  constrained  by  tide  information  through  assimilation.  The  participating 
seven  models  are  listed  in  the  upper  part  of  Table  1.  As  can  be  seen  from  the  table,  nominal  model  resolu¬ 
tions  vary  from  7.5'  to  30'.  The  first  four  models  described  below  are  semiempirical  ("semi"  because  some 
of  their  small-scale  features  can  arise  from  their  underlying  prior  models),  while  the  last  three  are  barotropic 
hydrodynamical  models  constrained  by  assimilation.  All  assimilative  models  include  the  constituents  M2,  S2, 
N2,  K2,  O, ,  P, ,  Q1 ,  and  K, .  In  addition,  some  also  provide  M4,  MS4,  and  possibly  other  compound  tides. 

2.1.1.  GOT4.8 

The  GOT4.8  model  was  developed  at  Goddard  Space  Flight  Center  and  follows  a  long  series  of  similar  efforts 
starting  with  Schrama  and  Ray  [1994]  and  Ray  [1999],  Like  several  of  the  other  models,  it  is  the  result  of  an 
empirical  harmonic  analysis  of  satellite  altimetry  relative  to  an  adopted  prior  model.  In  this  case  the  prior 


STAMMER  ETAL. 


©2014.  American  Geophysical  Union.  All  Rights  Reserved. 


245 


IgAGU  Reviews  of  Geophysics 


10.1 002/201 4RG000450 


was  a  collection  of  both  global  and  regional  models  blended  at  mutual  boundaries.  Of  the  modern  mod¬ 
els  examined  here  GOT4.8  has  the  coarsest  spatial  resolution  of  0.5°,  and  it  is  therefore  expected  to  be  less 
accurate  in  near-coastal  water. 

In  the  deep  ocean  between  latitudes  ±66  °,  the  GOT4.8  solution  was  based  only  on  data  from  the  T/P  satel¬ 
lite,  including  its  extended  mission  along  an  interlaced  ground  track;  no  Jason  data  were  used,  in  part  to 
facilitate  some  consistency  checks  between  T/P  and  Jason  [Ray,  201 3,  Appendix  A].  In  shallow  seas  and  also 
in  deep  water  poleward  of  66°,  data  from  Geosat  Follow-On,  ERS-1,  and  ERS-2  were  used.  A  small  amount 
of  early  ICESat  data  were  used  in  the  Weddell  and  Ross  seas.  The  load  tide,  i.e.,  the  Earth's  crustal  deforma¬ 
tion  under  the  weight  of  the  overlying  ocean  tide,  was  accounted  for  by  an  iterative  method  [Cartwright 
and  Ray,  1991],  Note  that  the  earlier  version  GOT4.7  is  identical  to  GOT4.8  except  for  the  S2  constituent.  The 
later  version  accounts  for  a  small  error  caused  by  the  S2  atmospheric  tide  in  the  T/P  dry-troposphere  cor¬ 
rection.  Correcting  for  this  small  error  makes  a  noticeable  improvement  in  the  final  S2  ocean  tide  solution 
[Ray,  201 3], 

Barotropic  currents  were  derived  from  the  GOT4.8  elevations  through  a  least  squares  solution  of  the 
two-dimensional  momentum  and  continuity  equations  [Ray,  2001]. 

2.1.2.  OSU12 

OSU12  is  a  global  ocean  tide  model  developed  at  the  Ohio  State  University  (OSU)  [Fok,  2012].  It  is  an  empir¬ 
ical  ocean  tide  model  determined  using  an  extended  orthogonal  representation  of  the  response  method 
[Groves  and  Reynolds,  1 975].  A  spatiotemporal  weighting  algorithm  was  employed  for  the  solution  to  com¬ 
bine  multimission  radar  altimeter  data  with  distinct  tidal  aliasing  characteristics,  spatial  and  temporal 
coverage,  and  accuracy.  These  data  included  T/P,  Jason-1,  Geosat  Follow-On  (GFO),  and  Envisat,  covering 
the  period  October  1 992  through  January  2009,  with  latitude  coverage  limited  within  ±66°  and  load  tide 
corrections  based  on  the  NA099.b  model  of  Matsumoto  et  al.  [2000].  Outside  the  ±66°  latitude  bounds, 
the  current  OSU12v1 .0  model  is  patched  using  the  values  from  the  GOT4.7  model;  however,  efforts  are 
underway  to  extend  the  OSU1 2  tide  model  to  cover  the  polar  oceans.  The  OSU1 2  model  is  released  with 
eight  major  semidiurnal  and  diurnal  tides,  together  with  the  St  and  M4  tides.  The  model  is  available  at 
http://geodeticscience.Org/oceantides/OSU12v1.0,  where  future  updates  to  the  model  will  also  be  hosted. 

2.1.3.  DTU10 

The  DTU10  global  ocean  tide  model  [Cheng  and  Andersen,  2011]  is  based  on  an  empirical  correction  to  the 
global  tide  model  FES2004  [Lyard  et  al.,  2006],  in  which  the  largest  residual  tides  were  determined  using  the 
response  method.  In  order  to  improve  shallow-water  tides,  4  years  of  theTopex/Jason-1  interleaved  mission 
data  were  introduced.  Combined  Envisat,  GEOSAT  Follow-On,  and  ERS-2  data  sets  were  introduced  outside 
the  ±66°  parallels  to  derive  tides  in  the  polar  seas.  As  Envisat  and  ERS-2  were  Sun-synchronous  satellites, 
diurnal  tides  could  be  significantly  improved  by  removing  annual  sea  level  variations  prior  to  estimating 
the  residual  tides.  The  four  major  semidiurnal  and  diurnal  tides  were  interpolated  onto  a  regular  0.1 25°  by 
0.125°  grid  using  a  depth-dependent  interpolation  method  [Andersen,  1995, 1999].  The  load  model  used  is 
that  of  FES2004.  The  model  can  be  downloaded  via  ftp://ftp.space.dtu.dk/pub/DTU10/. 

2.1.4.  EOTIIa 

EOT1  la  [Savcenko  and  Bosch,  201 2]  is  the  latest  version  of  a  series  of  global  ocean  tide  models,  derived  at 
Deutsches  Geodatisches  Forschungsinstitut  (DGFI)  by  means  of  empirical  analysis  of  multimission  satel¬ 
lite  altimetry  data.  EOT1 1  is  used  for  reprocessing  of  GRACE  gravity  data  [Bosch  et  al.,  2009];  it  is  based  on 
sea  surface  heights  observed  by  Topex,  Jason-1,  Jason-2,  ERS-2,  and  Envisat  over  the  ice-free  oceans  cov¬ 
ering  altogether  the  period  October  1 992  to  March  2010.  Tidal  constituents  were  estimated  based  on  a 
residual  least  squares  harmonic  analysis  with  FES2004  taken  as  reference  model.  EOT1  la  is  provided  on 
the  same  regular  1  /8  °  grid  as  FES2004.  At  latitudes  poleward  of  ±81 .5  °  (the  latitude  limits  of  ERS-2  and 
Envisat)  EOT1 1  a  defaults  back  to  FES2004.  The  EOT1 1  a  methodology  accounted  for  correlations  among  the 
along-track  altimeter  observations,  applied  a  careful  editing  procedure,  and  used  a  variance  component 
estimate  in  order  to  achieve  a  relative  weighting  of  data  from  different  missions.  In  addition,  EOTIIa  solved 
for  mission-specific  offsets  to  account  for  residual  altimeter  biases.  An  analysis  of  the  variance-covariance 
matrix  indicated  only  small  correlations  among  all  constituents  commonly  estimated.  Therefore,  EOT1  la 
is  composed  of  individual  constituents  taken  from  different  solutions  with  varying  smoothing  scales,  with 
preference  given  to  those  solutions  leading  to  minimal  root-mean-square  (RMS)  differences  to  validation 
data.  The  load  tides  were  computed  following  the  algorithm  of  Cartwright  and  Ray  [1991].  Model  fields  can 
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be  downloaded  at  ftp://ftp.dgfi.badw.de/pub/EOT1 1  a/data  and  are  described  in  detail  by  Savcenko  and 
Bosch  [2012], 

2.1.5.  HAMTIDE 

The  HAMTIDE12  (HAM12,  for  short)  model  is  a  global  ocean  barotropic  data-assimilative  model  developed 
at  the  University  of  Hamburg  [Zahel,  1995;  Taguchi  et  al.,  201 4].  The  model  is  based  on  the  linearized  tidal 
hydrodynamic  equations  and  on  a  simple  harmonic  time  dependence  of  the  variables.  The  model  includes 
the  full  loading  and  self-attraction  (LSA)  effect  (also  known  as  SAL  [ Hendershott ,  1972]).  Finite  difference 
methods  were  adopted  for  numerically  solving  the  equations.  The  continuous  space-time  domain  was  con¬ 
verted  into  the  FES2004  [Lyard  et  al.,  2006]  discrete  grid  domain,  where  the  grid  spacing  becomes  smaller 
in  higher  latitudes.  The  bottom  topography  is  based  on  the  GEBCO  1 '  bathymetry  field.  EOT1 1  a  output  was 
used  between  74°N  and  84°S  as  a  constraint.  A  direct  optimization  method  was  used  for  finding  the  mini¬ 
mum  of  a  cost  function  containing  dynamical  residual  smoothing  terms,  which  are  the  first  and  the  second 
derivatives  of  the  model  residuals.  Thus,  the  model  does  not  smooth  control  variables  directly.  The  adopted 
numerical  assimilation  scheme  allows  some  compensation  for  model  deficiencies  [Taguchi  et  al.,  2014]. 

2.1.6.  FES2012 

FES2012  (FES12,  for  short)  is  the  most  recent  of  a  series  of  global  finite  element  tidal  solutions,  initiated  by  Le 
Provost  et  al.  [1 994].  The  FES1 2  atlas  [Carrere  et  al.,  2012]  is  built  upon  altimetry-derived  harmonic  constant 
assimilation  using  the  ensemble,  frequency  domain  SpEnOI  (Spectral  Ensemble  Optimal  Interpolation)  data 
assimilation  software  in  T-UGOm  (Toulouse-Unstructured  Grid  Ocean  model)  hydrodynamic  tidal  solutions. 
T-UGOm  is  a  2-D/3-D  ocean  hydrodynamics  model  based  on  unstructured  meshes.  It  has  been  carried  out  in 
time-stepping  mode  and  frequency  domain  mode.  Frequency  domain  2-D  governing  equations  are  derived 
from  the  classical  shallow-water  continuity  and  momentum  equations,  as  described  in  Lyard  et  al.  [2006]. 
The  horizontal  discretization  used  for  FES1 2  is  continuous  Lagrange  polynomial  second-order  interpolation 
(LGP2)  for  elevation  and  discontinuous  nonconforming  PI  (linear  approximation  with  nodes  located  at  ele¬ 
ment  side  midpoints)  for  tidal  currents.  Resolution  varies  from  a  few  kilometers  in  coastal  areas  up  to  about 
25  km  in  the  deep  ocean.  The  unstructured  grid  covers  the  global  ocean,  with  the  required  topography 
based  on  a  number  of  sources,  including  approximately  20  regional  terrain  models,  with  considerable  effort 
devoted  to  improving  shelf  and  ice-shelf  seas.  Ensembles  were  built  by  perturbing  the  main  tidal  parame¬ 
ters,  e.g.,  bathymetry,  bottom  friction,  and  internal  tide  drag.  The  density  of  data  used  in  assimilation  was 
tuned  as  a  function  of  depth.  In  addition  to  T/P,  Jason-1,  and  Jason-2  data  which  are  the  FES1 2  main  assim¬ 
ilation  data  sources,  interleaved  mission  data  have  been  added  in  some  regions  showing  small-scale  tidal 
patterns.  Similarly,  ERS-1,  ERS-2,  and  Envisat  data  have  been  used  in  the  Arctic.  GLORYS2-V1  (Global  Ocean 
Reanalysis  2  Version  1 )  reanalysis  [Ferry  et  al.,  2012]  was  used  to  remove  nontidal  annual  and  semiannual 
contributions  from  altimetry  harmonic  analysis,  and  GOT4.8  was  used  for  load  tide  corrections. 

2.1.7.  TPX08 

TPX08  is  the  most  recent  in  a  series  of  tidal  solutions  (elevations  plus  barotropic  currents)  produced  using 
the  representer-based  variational  scheme  described  by  Egbert  et  al.  [1 994]  and  Egbert  and  Erofeeva  [2002] 
to  assimilate  altimetry  (and  other)  data  into  a  global  shallow-water  model.  The  TPX08-atlas  solution  dis¬ 
cussed  in  this  paper  represents  a  further  refinement,  with  a  series  of  over  30  regional  assimilation  solutions 
incorporated  to  increase  resolution  in  coastal  areas  and  shallow  seas.  The  base  global  solution  has  a  resolu¬ 
tion  of  1/6°,  with  bathymetry  for  the  dynamical  model  derived  from  the  GEBCO  V  database  [GEBCO  Digital 
Atlas,  2003]  with  some  adjustments  around  Antarctica  south  of  57°S.  Primary  data  were  harmonically  ana¬ 
lyzed  along-track  tidal  constants  from  685  T/P  and  Jason  cycles  (up  to  2011),  with  load  tide  corrections 
based  on  the  earlier  assimilation  solution  TPX06.2.  In  the  deep  ocean,  harmonic  constants  approximately 
every  1 2  km  along  track  were  assimilated.  For  shallower  depths  harmonic  constants  for  every  along-track 
point  were  used,  along  with  data  from  1 14  cycles  of  T/P  on  its  interleaved  ground  track.  At  high  latitudes, 
tide  gauge  data  were  also  assimilated,  including  83  around  Antarctica  [King  and  Padman,  2005]  and  289  in 
the  Arctic  [Kowalik  and  Proshutinsky,  1 994;  G.  Kivman,  personal  communication,  2002].  For  the  nonlinear 
quarter-diurnal  constituents,  the  two-stage  scheme  of  Egbert  et  al.  [2010]  was  used.  Regional  solutions  were 
obtained  with  a  resolution  of  1/30°  for  33  rectangular  areas,  including  all  major  enclosed  or  semienclosed 
seas  and  most  coastal  areas  with  significant  continental  shelf  width.  Higher-resolution  local  bathymetry 
was  used  wherever  possible,  and  all  available  T/P-Jason  and  interleaved  data  were  assimilated.  At  high 
latitudes,  or  in  shallow  areas  with  sparse  T/P  coverage,  ERS/Envisat  data  were  also  used  for  lunar  tides. 
Additional  coastal  tide  gauges  were  used  in  some  regional  solutions  to  improve  estimates  in  local  bays 
(see  http://volkov.oce.orst.edu/tides/region.html).  The  regional  solutions  were  patched  into  the  global 
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TPX08  base  solution,  keeping  higher  resolution  in  coastal/shallow  areas  and  using  a  weighted  average 
of  regional  and  global  solutions  over  a  narrow  strip  for  a  smooth  transition  to  the  global  model  in  the 
open  ocean.  TPX08-atlas  is  available  in  a  multiresolution  version,  but  the  fixed  2'  grid  was  used  in  all  tests 
performed  here. 

2.2.  Historical  Tide  Models 

To  help  gauge  the  progress  that  has  occurred  over  the  years,  we  give  statistics  for  two  historical  models  in  a 
few  sections  below.  We  use  one  model  from  the  era  before  T/P  and  one  which  was  developed  shortly  after 
the  launch  of  T/P  and  which  was  recommended  by  Shum97.  Respective  models  are  listed  in  the  middle  part 
of  Table  1 . 

2.2.1.  NSWC 

The  most  widely  adopted  global  tidal  model  in  the  pre-T/P  era  was  developed  at  the  U.S.  Naval  Surface 
Weapons  Center  (NSWC)  by  Schwiderski  [1 979].  He  employed  a  large  database  of  coastal  and  island 
harmonic  constants  and  assimilated  them  into  a  1  °  barotropic  tidal  model  using  a  method  he  called  "hydro- 
dynamic  interpolation,"  which  seems  to  have  been  a  variation  of  what  is  now  understood  as  a  relaxation  or 
"nudging"  method  of  assimilation  [ Egbert  and  Bennett,  1 996].  Where  he  had  good  data  to  guide  the  model, 
his  results  were  quite  good.  However,  in  some  wide  expanses  of  the  open  ocean,  such  as  the  mid-Atlantic, 
errors  in  excess  of  1 0  cm  were  later  found  [e.g.,  Cartwright  and  Ray,  1 990]. 

2.2.2.  CSR3.0 

One  of  the  objectives  of  Shum97  was  to  select  a  tidal  model  for  subsequent  reprocessing  of  satellite 
altimeter  data.  Based  on  various  rankings,  they  recommended  the  CSR3.0  model,  developed  by  Eanes  and 
Bettadpur  [1996].  Those  authors  had  used  a  response  analysis  applied  to  2.4  years  of  T/P  data.  They  adopted 
a  prior  hydrodynamic  model  based  primarily  on  FES94  of  Le  Provost  et  at.  [1 994].  The  adjustment  to  FES94 
was  heavily  smoothed,  which  worked  well  over  the  open  ocean  but  was  probably  less  suitable  for  regions 
near  land,  a  point  that  becomes  clear  in  some  of  our  tests  below  (e.g.,  see  section  4).  But  in  the  deep  open 
ocean,  which  has  been  the  primary  domain  of  satellite  altimetry,  CSR3.0  represented  the  new  paradigm  of 
tidal  accuracy  made  possible  by  T/P. 

2.3.  Purely  Hydrodynamic  Models 

Besides  data-constrained  models,  purely  hydrodynamic  (forward)  tide  models  are  also  tested  in  this  paper 
with  respect  to  their  skill  in  simulating  the  M2  barotropic  tidal  currents  and  sea  surface  elevations.  All 
respective  models  are  listed  in  the  lower  part  of  Table  1.  The  first  four  of  the  participating  forward  models 
described  below  are  barotropic,  while  the  latter  two  are  baroclinic  ocean  general  circulation  models  with 
concurrent  atmospheric  and  tidal  forcing. 

2.3.1.  HIM 

The  barotropic  version  of  the  Hallberg  Isopycnal  Model  (HIM)  [Hallberg  and  Rhines,  1996]  has  been  uti¬ 
lized  in  a  number  of  tidal  studies,  beginning  with  Arbic  et  al.  [2004].  The  simulations  were  run  on  a 
latitude-longitude  grid  from  86°S  to  82°N,  with  an  artificial  wall  at  the  northern  boundary.  The  source 
topography  files  are  described  in  detail  in  Arbic  et  al.  [2004];  in  brief,  the  topography  equatorward  of  72°  is 
taken  from  Smith  and  Sandwell  [1 997],  with  specialized  Arctic  and  Antarctic  topographic  data  sets  blended 
in  at  higher  latitudes.  Similar  to  other  recent  forward  tide  models  [e.g.,  Jayne  and  St.  Laurent,  2001],  the  HIM 
simulations  employ  a  parameterized  topographic  wave  drag  to  represent  the  tidal  energy  loss  occurring  in 
regions  of  strong  internal  wave  generation  over  rough  topography  [Egbert  and  Ray,  2000].  As  in  Egbert  et 
al.  [2004],  the  HIM  simulations  utilize  an  iterative  method  to  compute  the  SAL  term.  The  particular  simula¬ 
tion  used  here  is  the  1/8°  simulation  of  Arbic  et  al.  [2008],  which  yields  more  accurate  tides  than  the  coarser 
resolution  runs  of  Arbic  et  at.  [2004]. 

2.3.2.  OTIS-ERB 

The  OTIS-ERB  solution  was  derived  by  time-stepping  the  nonlinear  shallow-water  equations,  with  forcing 
at  the  M2  and  K,  frequencies,  on  a  1/12°  finite  difference  grid  running  from  86°S  to  82°N.  The  numerical 
discretization  and  much  of  the  software  are  identical  to  that  used  for  the  data-assimilating  solution  TPX08. 
Dissipation  was  parameterized  as  a  sum  of  a  quadratic  term,  representing  bottom  friction  in  shallow  seas, 
and  a  linear  internal  wave  drag,  comparable  to  that  used  by  Jayne  and  St.  Laurent  [2001  ].  SAL  was  rigorously 
modeled  with  an  iterative  scheme,  first  suggested  by  Accad  and  Pekeris  [1978];  these  second-order  forcing 
terms  were  computed  by  high-degree  spherical  harmonic  series.  Further  details  on  the  implementation  and 
validation  of  the  OTIS-ERB  solution  are  provided  in  Egbert  etaL  [2004]. 
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2.3.3.  OTIS-GN 

The  OTIS-GN  model  presented  by  Green  and  Nycander  [201 3]  is  basically  the  same  model  as  OTIS-ERB  [Egbert 
et  at.,  2004]  with  a  resolution  of  1  /8°  but  with  a  few  additional  modifications.  Firstly,  the  domain  for  OTIS-GN 
spans  80°S  to  80°N,  with  elevations  from  TPXO.6.2  being  used  as  boundary  conditions.  The  second  major 
change  was  that  a  modified  tidal  conversion  parameterization  based  on  the  scheme  presented  by  Nycander 
[2005]  was  used.  The  original  parameterization  was  modified  to  describe  the  tidal  conversion  in  terms  of  a 
2x2  tensor  rather  than  a  two-dimensional  vector  (see  Green  and  Nycander  [201 3]  for  further  details). 

2.3.4.  STM-1  B 

STM-1  B  is  a  depth-averaged  hydrodynamic  tidal  model,  formulated  in  terms  of  the  barotropic  shallow-water 
equations  in  volume  transport  form.  A  spectral  approach  is  used,  in  which  the  equations  of  motion  are 
projected  onto  a  set  of  discrete  tidal  frequencies,  so  that  the  time  dependence  is  removed,  leaving  a  set 
of  partial  differential  equations  in  latitude  and  longitude.  Nonlinearity  in  the  momentum  equation  is 
neglected,  apart  from  the  parameterized  bottom  drag  term,  which  takes  the  standard  form  pcd|u|u.  An  iter¬ 
ative  scheme  is  required  to  account  for  the  parameterized  nonlinear  bottom  drag  and  to  allow  for  a  full 
implementation  of  the  SAL  term  using  spherical  harmonic  transforms  (which  correspond  to  a  dense  matrix 
operation).  A  frequency-dependent  linear  internal  tide  drag  parameterization  is  used  [Griffiths  and  Peltier, 
2009].  The  model  uses  a  rotated  pole  (passing  through  Greenland  and  Antarctica)  and  is  thus  truly  global. 
Further  details  of  this  iterative  numerical  scheme  are  given  by  Hilt  etal.  [2011].  Model  results  were  calculated 
with  a  grid  spacing  of  1/12°  (2160  x  4320  grid  points),  with  forcing  from  three  semidiurnal  constituents 
(M2,  S2,  and  N2)  and  two  diurnal  constituents  (K,  and  O,).  In  the  ocean  interior,  energy-conserving 
second-order  finite  differencing  is  used  on  an  Arakawa  C  grid.  Flowever,  the  coastal  boundary  conditions 
were  implemented  using  a  standard  C  grid  staircase  coastal  methodology,  which  degrades  the  overall 
accuracy  of  the  solution  to  first  order  in  grid  spacing  [Griffiths,  2013]. 

2.3.5.  STORMTIDE 

The  STORMTIDE  model  [Mutter  et  at.,  2012]  is  a  high-resolution  ocean  circulation  and  tide  model  which 
resolves  mesoscale  eddies  and  low-mode  internal  tides.  It  has  been  developed  in  the  framework  of  the 
German  consortium  project  STORM,  which  aims  for  a  climate  model  simulation  at  the  highest  resolution 
currently  possible  [von  Storch  et  at.,  2012].  The  ocean  model  includes  an  embedded  thermodynamic  sea  ice 
model  and  is  based  on  the  Max-Planck  Institute  Ocean  Model  [Marsland  et  at.,  2003]  but  formulated  on  a 
tripolar  grid.  Ocean  tides  are  forced  in  real  time  by  a  lunisolar  tidal  potential  of  second  degree,  described 
by  analytical  ephemerides.  Further,  the  model  is  forced  by  a  climatological  wind  forcing,  which  includes 
changes  on  a  daily  basis  with  a  365  day  cycle.  Sea  surface  temperature  and  salinity  values  are  restored  to 
monthly  climatological  values,  in  order  to  simulate  a  realistic  seasonal  cycle  of  mixed  layer  depths.  The  grid 
of  the  model  has  40  vertical  z  levels  and  a  horizontal  resolution  of  about  0.1  “.The  SAL  effect  is  considered  in 
parameterized  form.  The  model  includes  no  internal  wave  drag;  instead,  in  the  deep  ocean  the  four  largest 
tidal  constituents  (M2,  S2,  K,,  and  O, )  convert  1 .1  TW  of  tidal  energy  from  barotropic  to  baroclinic  tides 
[Muller,  201 3].  The  model  topography  is  based  on  the  bathymetry  SRTM30_plus  (v6)  [Becker  et  at.,  2009]. 

2.3.6.  HYCOM 

This  baroclinic  tide  model  employs  the  forward  (nonassimilative)  simulations  of  the  FHYbrid  Coordinate 
Ocean  Model  (HYCOM),  run  at  high  vertical  and  horizontal  resolution  and  forced  by  the  astronomical  tidal 
potential  as  well  as  by  wind  and  buoyancy  forcing.  Chassignet  et  al.  [2007],  Metzger  et  at.  [201 0],  and  refer¬ 
ences  therein  provide  detailed  descriptions  of  HYCOM,  while  Arbic  etal.  [2010,  2012],  respectively,  provide 
details  and  an  overview  of  the  implementation  of  tides  in  HYCOM.  The  HYCOM  simulations  used  here  have 
32  layers  in  the  vertical  direction  and  an  equatorial  horizontal  resolution  of  1/12.5°.  The  land-sea  bound¬ 
ary  is  set  at  the  1 0  m  isobath.  The  bottom  topography  is  derived  from  the  Naval  Research  Laboratory  Digital 
Bathymetry  Data  Base  (DBDB2),  which  has  a  resolution  of  2  min  and  is  available  online  at  http://www7320. 
nrlssc.navy.mil/DBDB2_WWW/.  Numerous  hand  edits  were  performed  to  improve  coastlines  and  sill  depths 
in  key  straits  and  passages.  The  model  utilizes  a  "tripolar"  grid  and  is  therefore  truly  global.  The  SAL  effect 
is  currently  treated  with  the  simplified  scalar  approximation  of  Accad  and  Pekeris  [1978].  As  in  the  forward 
barotropic  studies  described  earlier,  a  parameterized  topographic  wave  drag  is  employed  in  the  HYCOM 
simulations  to  account  for  the  breaking  of  unresolved  high  vertical  mode  internal  waves  (see  Arbic  et  al. 
[2010]  for  more  discussion).  The  simulations  were  run  interannually  over  the  period  July  2003  to  December 
2010  using  6-hourly  Fleet  Numerical  Oceanography  Center  Navy  Operational  Global  Atmospheric  Prediction 
System  [Rosmond  et  at.,  2002]  atmospheric  forcing  with  wind  speeds  scaled  to  be  consistent  with  QuikSCAT 
[e.g.,  Liu  and  Xie,  2006]  observations.  The  HYCOM  surface  tidal  elevations  used  here  were  taken  from  Shriver 
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Table  2.  Standard  Deviations  of  Differences  Between  Tidal  Models  and  Their  RSSa 


Water  Depth 

Q, 

o. 

Pi 

Ki 

n2 

m2 

s2 

k2 

RSS 

Global 

0.22 

0.60 

0.32 

1.02 

0.50 

1.84 

0.86 

0.34 

2.458 

>  1000  m 

0.15 

0.23 

0.14 

0.30 

0.22 

0.30 

0.24 

0.19 

0.646 

<  1000  m 

0.23 

0.66 

0.38 

1.12 

0.78 

3.39 

1.50 

0.50 

4.061 

aOnly  data-constrained  models  listed  in  Table  1  are  included  in  this  comparison. 
RSS,  root-sum-square. 


etal.  [201 2]  and  are  based  on  an  analysis  of  year  2006.  For  the  barotropic  tidal  current  analysis  in  section  8, 
just  1  month  of  model  output  (September  2004,  output  at  hourly  intervals)  was  used  owing  to  the  enormous 
amount  of  storage  needed  for  a  high-resolution  three-dimensional  global  model. 


3.  Elevation  Differences  Between  Modern  Tide  Models 


Following  Andersen  etal.  [1995],  the  standard  deviation  (SD)  of  the  seven  tide  models  listed  in  the  top  part 
of  Table  1  was  computed  with  respect  to  elevations  r],  =  Ce~lat,  where  C  is  a  time-independent  amplitude 
of  a  tide  (e.g.,  M2  tide)  at  a  wet  grid  point  j,  a  denotes  its  frequency,  and  /  =  y—l.  As  a  first  step  of  the 
computation,  the  mean  elevation  of  each  tidal  constituent  across  all  models  ( N  =  7)  was  computed  at  every 
grid  point  according  to 

N 

4mean  ~  71  —  "mean  (COS  Gmean  +  /Sin  Gmean)  C  -  (1) 

;=i 

Flere  C  =  Hj  (cos  Gj  +  / sin  Gj)  with  amplitude  Hj  and  Greenwich  phase  lag  Gy  (hereafter,  "Greenwich"  is 
dropped  and  the  term  "phase  lag"  is  understood  to  be  in  reference  to  Greenwich,  as  is  standard  in  tidal 
studies).  The  SD  between  all  involved  models  can  then  be  computed  for  each  constituent  according  to 

/l  n  i  T  2 

-^tide  —  (  1u  T  /  —  *7mean))  ^ 

\  v  n= 1  •'O 

N  \  V2 

i  Z  T  [("->  C0S(Gn)  -  "mean  COS(Gmean))2  +  ("„  Sin(G„)  -  Hmean  Sin(Gmean))2]  )  ,  (2) 

n=1  / 


where  Hn  and  Gn  are  the  amplitude  and  Greenwich  phase  lag  of  a  constituent  given  by  each  model,  respec¬ 
tively,  and  Hmean  and  Gmean  are  the  amplitude  and  Greenwich  phase  lag  for  the  averages  of  equation  (1 ).  The 
SD  of  all  present  models  for  the  tide  (M2,  for  example)  is  computed  according  to 


SDtide  =  (  JX(SD‘deCOS(^)): 


k=  1 


1/2 


(3) 


Flere  k  is  the  grid  point  index,  P  is  the  number  of  all  elevation  grid  (wet)  points,  and  (pk  is  the  latitude  of  grid 
point  k. 

Computations  of  SDtide  were  performed  for  the  tidal  constituents  M2,  S2,  N2,  K2,  K,,  O,,  P,,  and  Q1  after 
regridding  bilinearly  the  fields  of  every  model  to  a  common  0.5°  grid.  Results  are  summarized  in  Table  2 
for  the  global  ocean  and  separately  for  the  deep  ocean  (>  1 000  m  depth)  and  the  shallow  seas  (<  1 000  m 
depth),  respectively.  As  an  example,  the  SD  of  all  models  is  1 .87  cm  for  M2  which  is  1 3%  lower  than  the 
2.29  cm  obtained  by  Shum97.  In  the  case  of  K,,  the  SD  of  1 .02  cm  exceeds  the  value  of  0.89  cm  obtained  by 
Shum97,  suggesting  that  the  present  models  diverge  more  from  each  other  than  was  the  case  previously. 
Flowever,  Shum97  did  not  include  polar  regions  and  was  significantly  less  inclusive  in  shallow  seas,  both 
areas  where  quite  large  SD  values  (often  >  2  cm  for  both  M2  and  1^  tides)  occur.  Moreover,  present  mod¬ 
els  agree  much  better  with  each  other  over  the  deep  ocean  (SD  =  0.3  cm)  as  compared  to  Shum97,  who 
reported  a  SD  of  0.5-1 .0  cm  for  both  M2  and  K,  tides  in  the  deep  sea  (>  1 000  m). 

Spatial  patterns  of  the  SD  of  model  differences  are  displayed  in  Figures  1  (top)  and  1  (bottom)  for  M2 
and  K, ,  respectively.  Clearly,  the  SD  of  both  constituents  is  large  in  polar  regions,  especially  in  the  Ross  and 
Weddell  Seas  in  the  Antarctic  and  in  Baffin  Bay,  between  Greenland  and  Ellesmere  Island.  The  same  holds 
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Figure  1.  Standard  deviations  for  (top)  M2  and  (bottom)  K,  from  the  seven  modern  data-constrained  models  listed 
in  Table  1. 


for  the  other  tidal  constituents  (not  shown)  and  presumably  results  from  the  lack  of  observations  there. 
However,  we  note  that  the  largest  discrepancies  between  models  can  be  found  in  the  coastal  Arctic  regions 
for  which  most  of  the  tide  gauge  information  is  available,  not  in  the  center  of  the  Arctic  where  tidal  ampli¬ 
tudes  are  much  smaller  (see  also  section  5.1).  In  the  Antarctic  the  largest  model  discrepancies  are  located 
under  the  ice  shelves;  an  exception  is  K,  for  which  the  largest  errors  are  in  the  open  continental  shelf  areas 
of  the  western  Ross  and  western  Weddell  Seas.  The  SD  of  semidiurnal  tides  are  also  enhanced  in  the  west¬ 
ern  portions  of  the  Pacific  and  Indian  Oceans,  characterized  by  complex  topography  and  in  some  cases  also 
extreme  tidal  amplitudes  (e.g.,  along  the  northern  Australian  continent,  in  the  Mozambique  Channel,  and  in 
the  Atlantic  on  the  Patagonian  Shelf),  along  the  northeastern  coast  of  Brazil  and  along  the  track  of  the  Gulf 
Stream,  where  the  magnitude  of  the  SD  is  between  0.3  and  1.2  cm  for  M2  and  S2  tides  and  0.3  and  0.6  cm 
for  N2  and  K2  tides.  Large  SD  of  diurnal  tides  resides  in  the  most  northern  part  of  the  Pacific  Ocean  with 
magnitudes  of  0.3  to  1 .2  cm  for  K,  and  O,  tides  and  0.3  to  0.6  cm  for  P1  and  Q,  tides.  Large  SD  is  also  found 
in  the  western  parts  of  the  Pacific  Ocean  and  in  the  Arabian  Sea.  In  the  Atlantic  Ocean,  SD  for  diurnal  tides 
is  small  everywhere  owing  to  the  generally  small  amplitudes  of  diurnal  tides.  Enhanced  model  discrepancies 
in  eddy-rich  western  boundary  current  regions  presumably  point  toward  eddy  energy  aliased  into  tide  mod¬ 
els  there.  We  also  note  that  intermodel  discrepancies  appear  to  be  substantially  elevated  in  shelf  regions 
relative  to  the  open  ocean  where  an  overall  RSS  value  over  the  eight  constituents  is  less  than  1  cm. 

4.  Tests  of  Modern  Tide  Models  Against  Tide  Gauge  Data 

Comparisons  with  in  situ  "ground  truth"  tidal  measurements  form  one  of  the  most  straightforward,  yet  pow¬ 
erful,  tests  of  model  output.  The  test  data  used  here  are  considerably  improved  and  extended  over  those 
used  by  Shum97.  Three  primary  data  sets  are  used,  one  based  in  deep  water,  one  in  shelf  water  and  shal¬ 
low  seas,  and  one  along  coastlines  of  continents.  (Polar  regions  are  reserved  for  section  5.)  The  philosophy 
for  separating  the  tests  in  this  way  rests  with  the  different  role  that  satellite  altimetry  plays  in  each  regime. 
Deep-ocean  tides,  for  example,  are  determined  primarily  by  direct  constraints  from  satellite  altimeter  mea¬ 
surements,  especially  from  Topex/Poseidon  and  the  Jason  follow-on  missions.  Similarly,  shelf  tides  can 
also  be  determined  in  large  part  by  altimetry,  especially  multimission  altimetry  [e.g.,  Cheng  and  Andersen, 
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201 1],  although  the  determination  is 
more  difficult  because  of  the  shorter 
scale  of  shelf  tides,  necessitating  hydro- 
dynamic  assimilation  (or  the  use  of 
hydrodynamic  priors  for  empirical 
models).  For  truly  coastal  tides,  how¬ 
ever,  altimetry  plays  little  role  except 
as  distant  constraints,  and  accuracies 
generally  must  rely  on  hydrodynamic 
modeling  and  an  accurate  bathymetry 
[e.g.,  Ray  et  a!.,  2011]. 

The  following  comparisons  test  indi¬ 
vidual  tidal  constituents,  focusing  on 
major  constituents  that  are  released 
by  all  or  most  models  under  study.  At 
each  comparison  station,  tidal  con¬ 
stants  from  the  models  are  interpolated 
to  the  station  location  by  bilinear 
interpolation.  Station  locations  are 
shown  in  Figure  2,  separated  into  (top) 
deep,  (middle)  shelf,  and  (bottom) 
coastal  stations. 


Testing  individual  tidal  constituents 
does  have  one  important  limitation. 
Many  users  of  tidal  models  require 
only  accurate  height  predictions,  and 
they  are  less  concerned  with  individ¬ 
ual  constituents.  Errors  in  a  few  major 
constituents  do  not  give  a  complete 
accounting  of  errors  in  height  pre¬ 
diction.  This  is  probably  not  a  great 
concern  for  deep  water,  since  forming 
an  RSS  across  the  eight  major  constituents  should  yield  a  rough  guide  to  prediction  error.  In  shallow  water, 
however,  where  dozens  (or  more)  compound  tides  may  be  required  for  accurate  prediction,  our  reliance 
on  a  few  constituents — and  indeed  the  tested  models'  similar  reliance  on  the  same — is  inadequate.  These 
matters  are  discussed  further  by  Ray  et  at.  [2011],  who  attempt  to  gauge  the  errors  of  omission  in  tidal 
prediction.  For  present  purposes  the  main  implication  is  that  our  statistics  in  sections  4.2  and  4.3  must  be 
considered  optimistic  as  a  guide  to  overall  tidal  height  prediction.  In  some  shallow-water  regions  it  is  pos¬ 
sible  that  the  omission  error  from  missing  constituents  is  comparable  to,  or  exceeds,  the  error  in  the  M2 
constituent,  which  normally  dominates  the  error  spectrum. 


Figure  2.  Station  locations  for  the  "ground  truth"  tidal  data  used  in 
section  4.  (top)  Deep-water  stations,  (middle)  Shelf  stations,  (bottom) 
Coastal  stations. 


4.1.  Deep-Ocean  Gauges 

The  deep-ocean  test  data  consist  of  harmonic  constants  at  151  sites,  all  deduced  from  long  time  series 
of  pressure  measurements  taken  on  the  seafloor  (Figure  2,  top).  All  are  independent  of  the  models  being 
tested.  Use  of  open-ocean  bottom  pressure  avoids  local  tide  effects  that  sometimes  can  occur  even  at 
small  islands,  so  the  data  are  thus  comparable  to  the  tidal  fields  measured  by  satellite.  Compilation  of  these 
data  is  described  by  Ray  [2013],  who  provides  details  regarding  station  selection  criteria,  analysis  meth¬ 
ods,  and  various  corrections.  Preference  was  given  to  long  time  series;  the  median  length  was  567  days.  The 
solar  semidiurnal  S2  tide  in  bottom  pressure  was  corrected  for  the  presence  of  the  atmospheric  barometric 
tide,  using  a  climatological  model  based  on  operational  surface  pressures  from  the  European  Centre  for 
Medium-Range  Weather  Forecasts  (ECMWF).  Conversion  of  bottom  pressure  amplitudes  to  equivalent  sea 
surface  heights  was  based  on  a  climatological  ocean  density  at  each  station  location. 

For  all  151  stations  the  signals  and  RMS  differences  for  major  constituents  are  summarized  in  Table  3,  where 
the  first  two  models  listed  are  the  two  historical  ones  dating  from  before  and  shortly  after  the  launch 
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Table  3.  Tide  Signal  and  RMS  Model  Differences  (cm)  Against  Deep-Ocean  Bottom  Pressure  Recorder  (BPR)  Stations3 


Ql 

o. 

Pi 

K, 

n2 

M2 

s2 

k2 

RSS 

M4 

Signal 

1.8 

8.8 

4.0 

12.5 

6.4 

30.2 

11.2 

3.1 

0.2 

Pre-Topex/Poseidon 

NSWC 

0.290 

0.874 

0.638 

1.292 

1.153 

4.268 

1.779 

0.660 

5.106 

Early  Topex 

CSR3.0 

0.230 

0.502 

0.252 

0.585 

0.375 

0.923 

0.607 

0.470 

1.514 

Modem  models 

GOT4.8 

0.165 

0.296 

0.234 

0.423 

0.252 

0.510 

0.369 

0.209 

0.923 

0.089 

OSU12 

0.304 

0.369 

0.194 

0.430 

0.441 

0.578 

0.940 

0.287 

1.394 

0.137 

DTU10 

0.226 

0.277 

0.292 

0.449 

0.274 

0.613 

0.415 

0.383 

1.088 

0.089 

EOTIIa 

0.232 

0.317 

0.224 

0.404 

0.335 

0.564 

0.428 

0.365 

1.056 

0.282 

HAM12 

0.160 

0.317 

0.199 

0.373 

0.245 

0.513 

0.397 

0.176 

0.904 

FES12 

0.216 

0.309 

0.355 

0.471 

0.342 

0.658 

0.407 

0.223 

1.120 

0.115 

TPX08 

0.153 

0.310 

0.181 

0.442 

0.201 

0.523 

0.338 

0.151 

0.893 

0.069 

Bootstrap  a 

0.013 

0.018 

0.012 

0.018 

0.018 

0.034 

0.024 

0.018 

0.008 

aModels  listed  in  order  of  increasing  spatial  resolution.  Bold  font  marks  models  within  1  cr  of  the  lowest  RMS  in 
any  given  constituent,  where  cr  is  taken  from  the  median  bootstrap  of  all  models  except  NSWC,  CSR,  and  TPX08. 


of  T/P  in  1 992.  The  signal  is  defined  by  y  ^xa:mn/2,  where  the  /\station  values  are  the  amplitudes  at  the 
respective  stations  and  the  overbar  represents  an  average  over  all  stations.  The  RMS  difference  is  given  by 
([^station  cos(rwf  -  </»station)  -  ^modei  cos(ruf  -  (/'model)]2)1  /2 -  where  here  the  overbar  is  computed  over  one  full 
cycle  of  the  constituent  in  question  (e.g.,  cot  varying  from  0  to  2k)  as  well  as  over  all  the  station  locations,  co  is 
frequency,  f  is  time,  tp  denotes  phase  lag,  and  the  "model"  subscript  denotes  model  values.  A  striking  feature 
of  the  table  is  the  very  low  RMS  for  all  modern  models:  half  a  centimeter  or  better  for  all  constituents,  repre¬ 
senting  relative  errors  (a  ratio  of  RMS  difference  to  signal)  ranging  from  about  2%  for  M2  up  to  10%  for  the 
smaller  K2,  P,,  and  constituents.  In  contrast,  the  comparable  deep-ocean  data  set  used  by  Shum97 — a 
102-station  set  compiled  primarily  by  Christian  Le  Provost  and  David  Cartwright — showed  much  larger 
RMS  differences,  e.g.,  1.6  cm  for  M2.  We  note,  however,  that  much  of  the  improvement  here  comes  from 
the  improved  test  data  set,  not  improvements  in  the  models.  For  example,  for  GOT4.8  the  old  102-station 
test  set  gives  an  RMS  of  1 .45  cm  [Ray,  2013]  versus  the  0.51  cm  RMS  value  here.  The  new  data  set  benefits 
from  many  newer  multiyear  bottom  pressure  time  series  now  available,  including  those  of  the  international 
tsunami  warning  network,  while  the  old  data  set  had  many  very  short  (30  days)  bottom  pressure  series  and 
some  island  data  that  were  evidently  biased  by  local  harbor  or  lagoon  perturbations.  Nonetheless,  Table  3 
also  shows  good  improvement  over  the  historical  models,  with  the  RMS  for  CSR3.0  M2  nearly  double  that 
of  the  modern  models,  and  NSWC  almost  an  order  of  magnitude  worse.  Most  of  the  other  constituents  also 
show  good  improvement  over  the  1 997-era  "best"  model. 

The  final  line  of  Table  3  is  our  attempt  to  gauge  statistical  uncertainty  in  the  computed  RMS  values  in  order 
to  assess  the  significance  of  very  small  differences  among  models.  For  six  of  the  models  we  computed  a 
bootstrap  estimate  of  the  standard  deviation  of  the  RMS  values  (resampling  with  replacement),  and  the 
quoted  entries  in  the  table  are  the  median  standard  deviations  for  each  constituent.  We  consider  any  mod¬ 
els  having  RMS  values  within  an  interval  less  than  this  bootstrap  a  to  be  equivalent,  and  bold  font  in  the 
table  marks  all  models  within  la  of  the  lowest  RMS  in  any  given  constituent. 

Four  of  the  modern  models  in  Table  3  have  at  least  one  bold  entry,  indicating  a  best  ranking  in  that  con¬ 
stituent.  TXPO.8  has  the  most  such  entries,  with  best  results  in  all  four  semidiurnal  constituents;  only  its 
K,  shows  somewhat  higher  RMS  in  comparison.  The  S2  constituent  for  OSU1 2  is  anomalous  with  RMS  of 
0.94  cm.  In  this  case  we  speculate  a  problem  with  an  altimeter  correction,  since  many  of  these,  including 
atmospheric  delay  corrections  and  atmospheric  loading  corrections,  can  have  significant  power  at  the  S2 
frequency.  An  overweighting  of  Sun-synchronous  altimeter  data,  which  can  corrupt  S2  estimates,  is  another 
potential  explanation  for  the  anomalous  OSU1 2  S2  RMS  value. 

An  RMS  difference  is  only  one  measure  of  discrepancy.  To  test  the  robustness  of  these  results,  we  have  also 
computed  the  median  absolute  difference  between  in-phase  and  quadrature  components  of  models  and 
test  stations.  Such  median  differences  should  be  less  sensitive  to  a  few  bad  comparisons,  although  there  are 
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very  few  of  those  in  the  deep-ocean  stations.  The  median  results  are  not  shown  here,  but  they  are  broadly 
consistent  with  results  of  Table  3.  In  particular,  DTU10  shows  the  best  ranking  in  O,  and  HAM1 2  by  far  the 
best  in  Kv  while  TPX08  is  again  best  in  all  four  semidiurnal  constituents  as  well  as  in  Q,  and  P1  (with  ties  if 
statistical  variability  is  allowed  for). 

4.2.  Shallow/Shelf  Sea  Gauges 

The  test  data  of  this  section  consist  of  harmonic  constants  from  195  stations  located  in  water  generally  shal¬ 
lower  than  200  m  but  not  along  continental  coastlines  (Figure  2,  middle).  These  stations  therefore  represent 
primarily  shelf  tides.  Given  generally  larger  tidal  amplitudes  and  shorter  wavelengths  in  these  regions  rel¬ 
ative  to  tides  in  the  deep  ocean,  we  must  anticipate  larger  discrepancies  with  the  in  situ  data  and,  as  seen 
in  section  3,  among  the  models  themselves.  Nonetheless,  shelf  tides  can  be  constrained  by  altimetric  data, 
even  if  spatial  resolution  is  challenging,  and  extending  satellite  altimetry  into  shallow-water  regions  is  a 
topic  of  growing  interest  [e.g.,  Vignudelli  eta!.,  201 1],  We  therefore  wished  to  design  a  test  appropriate  to 
shallow,  but  not  coastal,  water.  Thus,  in  this  section  we  avoid  tide  gauges  located  along  coastlines  of  large 
land  masses,  although  we  do  employ  gauges  at  small  offshore  islands;  this  is  especially  the  case  for  our  data 
near  the  Australian  continent,  which  are  all  offshore  island  gauges.  Otherwise,  the  great  majority  of  stations 
are  from  bottom  pressure  shelf  measurements. 

The  tide  gauge  data  here  build  on  the  data  set  compiled  and  used  by  Rayetal.  [2011]  but  now  augmented 
to  a  total  of  195  stations.  Newly  added  stations  are  extracted  from  recent  publications,  such  as  two  new 
stations  on  the  Patagonian  Shelf  off  Tierra  del  Fuego  [ Richter  etal.,  2012],  or  from  previously  overlooked 
publications,  such  as  data  from  the  Gulf  of  Gabes  [Sammari  et  al.,  2006]  and  South  China  Sea  [Beardsley  et 
al.,  2004],  as  well  as  from  unpublished  data  (e.g.,  data  from  the  Savu  Sea  provided  by  Dr.  Janet  Sprintall, 
October  2009).  Nevertheless,  the  spatial  coverage  of  the  data  is  still  far  from  adequate,  as  is  evident  from 
Figure  2  (middle). 

These  shelf  data  are  themselves  much  less  accurate  than  the  deep-ocean  data  of  the  previous  section — they 
include  shorter  time  series,  stations  less  well  documented,  and  stations  displaying  considerably  higher 
differences  with  models  but  not  obviously  identifiable  as  being  in  error.  In  a  few  isolated  cases  we  have 
used  anomalous  discrepancies  with  our  seven  primary  models  to  identify  and  sometimes  correct  problem 
stations.  One  case  involved  new  bottom  pressure  data  from  the  northern  Adriatic  Sea  [Book  et  al.,  2009];  cor¬ 
rected  harmonic  constants  were  kindly  provided  by  Dr.  Jeffrey  Book.  Another  case  involved  one  of  the  few 
stations  from  the  interior  of  the  Yellow  Sea,  with  tidal  constants  published  by  Teague  et  al.  [1 998];  their 
Station  C  is  located  in  the  middle  of  the  Yellow  Sea  and  should  be  well  observed  by  satellite,  but  the 
reported  phase  lag  for  M2  of  1 69.3°  is  in  poor  agreement  with  our  models  (with  values  1 82.7°,  1 82.0°,  1 81 .8°, 
183.7°,  181.8°,  182.3°,  and  182.2°).  The  cause  of  this  problem  has  not  been  resolved,  but  ultimately  we 
removed  this  station  from  our  test  set.  (It  is  worth  noting  that  Yao  et  al.  [2012],  in  their  regional  modeling  of 
the  Yellow  Sea  tides,  obtained  a  phase  of  1 81  °  in  their  unconstrained  prior  solution  but  1 72°  after  assimilat¬ 
ing  the  Station  C  data;  we  surmise  that  their  prior  may  have  been  the  more  accurate  at  that  location.)  There 
are  a  few  similar  cases.  In  all  cases  we  had  confidence  that  any  rejected  tide  gauge  station  was  in  a  location 
that  should  have  been  well  observed  by  altimetry,  and  yet  models  differed  unusually  from  gauge  estimates; 
otherwise  we  kept  the  gauge  even  if  differences  were  large. 

Of  the  195  test  stations,  76  are  on  the  northwest  European  Shelf.  To  avoid  overweighting  that  geographic 
region  in  our  global  statistics,  and  to  compare  consistency  of  statistics,  Table  4  gives  RMS  gauge  model  dif¬ 
ferences  split  between  the  European  Shelf  and  the  remainder  of  the  globe.  (Unlike  Table  3,  Table  4  could 
not  include  the  NSWC  model  because  that  model  was  undefined  for  the  Mediterranean  Sea  and  the  Strait 
of  Korea;  but  leaving  out  those  stations,  its  RMS  for  M2  is  33.1  cm,  far  larger  than  any  entry  in  the  table.)  Of 
the  195  stations,  10  of  them  were  assimilated  by  TPXO.8,  mostly  in  the  region  around  Australia.  None  of  the 
European  tide  gauges  were  assimilated,  so  the  assimilation  affects  only  the  lower  portion  of  Table  4.  None  of 
the  other  models  assimilated  any  of  our  test  stations. 

Unlike  the  deep-ocean  case,  now  all  seven  modern  models  have  at  least  one  constituent  marked  as  best  in 
Table  4.  As  before,  TPXO.8  has  the  most  such  entries,  which  may  partly  reflect  its  assimilation  of  1 0  of  the 
stations,  although  its  performance  is  also  excellent  on  the  European  Shelf  where  it  assimilated  no  stations. 
For  M2  on  the  European  Shelf,  the  models  cluster  into  two  groups  with  RMS  values  either  less  than  4  cm  or 
greater  than  5  cm.  The  best  models  are  DTU10,  PIAM12,  FES12,  and  TPX08.  Elsewhere,  the  best  M2  is  DTU10, 
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Table  4.  Tide  Signal  and  RMS  Model  Differences  (cm)  Against  Shelf  Water  Tide  Stations3 


Ql 

o, 

Pi 

K, 

n2 

M2 

s2 

k2 

RSS 

m4 

MS4 

European  Shelf 

Signal 

2.1 

5.6 

1.5 

5.5 

17.0 

87.7 

30.0 

7.3 

4.7 

2.9 

CSR3.0 

1.07 

2.10 

0.54 

1.67 

5.09 

13.34 

5.61 

2.38 

15.78 

GOT4.8 

0.93 

0.92 

0.55 

1.30 

1.97 

5.87 

2.51 

1.09 

7.04 

2.80 

OSU12 

1.11 

1.24 

0.69 

1.53 

1.77 

5.04 

4.04 

1.29 

7.22 

2.10 

DTU10 

0.83 

0.81 

0.51 

1.27 

2.17 

3.50 

2.38 

0.92 

5.17 

2.74 

EOTIIa 

0.85 

0.83 

0.50 

1.24 

2.13 

5.53 

3.43 

1.13 

7.17 

3.16 

HAM  12 

0.92 

1.96 

0.47 

1.14 

1.65 

3.11 

2.64 

0.92 

5.14 

FES12 

0.88 

0.82 

0.71 

1.19 

1.39 

3.71 

1.94 

0.63 

4.82 

2.22 

1.86 

TPX08 

0.88 

0.72 

0.46 

1.21 

1.58 

3.85 

1.70 

0.74 

4.87 

0.35 

1.16 

Bootstrap  a 

0.05 

0.08 

0.07 

0.11 

0.20 

0.70 

0.46 

0.14 

0.35 

Elsewhere 


Signal 

2.7 

11.0 

6.0 

17.3 

12.5 

53.6 

21.9 

7.1 

1.9 

1.5 

CSR3.0 

1.20 

2.44 

1.45 

4.96 

6.00 

23.11 

7.86 

2.65 

25.91 

GOT4.8 

0.69 

1.05 

0.92 

1.68 

1.97 

4.14 

2.93 

1.59 

6.11 

1.25 

OSU12 

1.00 

1.17 

0.84 

1.75 

1.87 

4.61 

3.00 

1.15 

6.42 

1.44 

DTU10 

0.82 

1.11 

1.06 

1.70 

1.80 

3.44 

2.39 

1.57 

5.40 

1.26 

EOTIIa 

0.73 

1.07 

0.78 

1.64 

1.86 

5.05 

3.39 

1.31 

6.87 

1.11 

HAM  12 

0.97 

1.16 

0.80 

2.02 

2.01 

3.89 

2.52 

1.44 

5.88 

FES12 

0.80 

1.00 

0.89 

1.51 

1.58 

3.33 

2.30 

1.02 

4.96 

0.98 

1.33 

TPX08b 

0.82 

1.00 

0.82 

1.47 

2.00 

3.50 

1.93 

1.12 

5.07 

0.88 

1.51 

Bootstrap  a 

0.08 

0.12 

0.08 

0.14 

0.17 

0.41 

0.21 

0.18 

0.12 

0.14 

aModels  listed  in  order  of  increasing  spatial  resolution.  Bold  font  marks  models  within  1  cr  of  the  low¬ 
est  RMS  in  any  given  constituent,  where  a  is  taken  from  the  median  bootstrap  of  all  models  except  CSR 
and  TPX08. 

bTPX08  assimilated  a  subset  of  the  test  stations. 


FES12,  HAM12,  and  TPXO.  GOT4.8  is  the  poorest  on  the  European  Shelf,  evidently  marred  by  its  coarse  0.5° 
resolution,  although  its  statistics  for  "Elsewhere"  are  relatively  good,  including  the  best  Q, . 

In  shallow  seas  it  is  unclear  whether  the  RMS  statistic  provides  a  meaningful  general  description  of  a  model. 
This  is  borne  out  by  Figure  3  which  shows  histograms  of  M2  vector  differences  (absolute  magnitudes) 
between  each  of  the  195  stations  and  two  models — one  of  the  best  (DTU10)  according  to  Table  4  and  one 
of  the  worst  (GOT4.8).  For  both  models,  most  station-model  differences  are  small,  less  than  2.5  cm,  but  there 
are  a  number  of  much  larger  differences  that  will  tend  to  skew,  rightly  or  wrongly,  an  RMS  statistic.  The  inter¬ 
pretation  is  that  in  most  shallow-water  locations  the  model  errors  in  M2  will  be  less  than  2.5  cm,  but  there 
are  still  some  locations  (not  always  evident  in  Figure  1)  where  errors  are  much  greater  (tens  of  centimeters). 
As  noted  above  the  error  can  also  be  assessed  by  the  median  absolute  differences  in  in-phase  and  quadra¬ 
ture  components.  Such  median  difference  results  for  M2  (not  shown)  vary  between  1 .1  cm  and  1 .5  cm,  with 
DTU10  and  FES12  slightly  superior.  We  emphasize  that  in  shallow  regions  one  single  statistic  cannot  ade¬ 
quately  capture  the  error  characteristics  of  these  models,  and  a  model  user  could  well  be  disappointed  by 
any  of  these  models  if  the  location  of  interest  is  one  of  the  outlier  stations  shown  in  Figure  3. 

4.3.  Coastal  Tide  Gauges 

Although  some  very  high-resolution  global  models  may  yield  adequate  tide  predictions  along  some  coast¬ 
lines,  they  cannot  be  expected  to  be  competitive  with  well-constructed  local  models  based  on  high-quality 
local  bathymetric  data  and  local  tidal  knowledge.  Nonetheless,  some  assessment  of  the  models  here  under 
consideration  may  be  of  value  to  those  wishing  to  rely  on  global  models.  For  example,  relying  on  global  tide 
models  is  not  uncommon  when  corrections  are  needed  for  satellite  altimetry,  since  local  models  are  often 
unavailable  while  the  global  models  are  readily  at  hand.  This  section  is  provided  primarily  to  stress  to  such 
users  that  tide  model  errors  can  be  potentially  very  large  in  the  near-coastal  zone. 

The  set  of  coastal  tide  gauges  used  in  this  section  is  a  subset  of  those  analyzed  by  Ponchautet  al.  [2001]. 
Those  authors  produced  high-quality  tidal  analyses  for  gauges  associated  with  the  World  Ocean  Circulation 
Experiment,  and  we  here  adopt  only  those  stations  situated  along  continental  coastlines  (or  nearby  islands, 
such  as  the  station  at  Zanzibar,  Tanzania).  After  removing  open-ocean  islands  from  the  Ponchaut  et  al.  data 
set  there  are  56  stations  left  (Figure  2,  bottom).  Although  this  number  is  small  compared  with  the  total 
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potential  data  set  of  coastal  stations, 
the  data  set  offers  advantages  of 
being  well  documented,  with  reliable 
tidal  harmonic  constants.  For  nearly 
all  tested  models,  some  station  loca¬ 
tions  were  deemed  by  the  model  to 
be  land  and  a  small  amount  of  extrap¬ 
olation  was  necessary  to  arrive  at  a 
valid  model  amplitude  and  phase.  Of 
the  56  stations,  two  of  them  (both  on 
the  Australian  coast)  were  assimilated 
by  TPXO.8. 

Table  5  shows  both  RMS  and  median 
absolute  differences  between  the  mod¬ 
els  and  the  56  coastal  stations.  In  this 
case,  even  more  so  than  with  the  shelf 
tides,  the  RMS  statistic  is  probably 
not  a  good  general  descriptor,  since 
one  or  two  poor  stations  can  domi¬ 
nate  the  results.  That  is  clearly  the  case 
for  the  M2  differences  of  FIAM1 2  and 
TPX08,  which  are  both  anomalously 
large.  TPX08  has  poor  agreement  at 
two  difficult  stations,  Zanzibar  and  Dar¬ 
win  (Australia),  while  HAM12  is  poor 
at  Zanzibar  and  Ketchikan  (Alaska). 
Another  station  which  has  poor  agree¬ 
ment  with  most  models  is  Duck  (North 
Carolina),  which  sits  on  a  barrier  island 
that  is  not  adequately  demarcated  in 
some  models,  resulting  in  errors,  some¬ 
times  on  both  sides  of  the  island — see 
Figure  4.  While  these  problem  stations 
skew  the  global  RMS  statistics,  they  are 
nonetheless  instructive  for  the  types  of  error  that  may  be  encountered  along  some  coastlines. 

The  more  robust  median  differences  of  Table  5  are  much  smaller  than  the  RMS  differences,  indicating  rea¬ 
sonably  accurate  model  values  for  the  majority  of  stations,  with  M2  component  differences  between  1.0  and 
1.5  cm  for  all  seven  modern  models  and  FI  AM  12  the  best  for  M2  notwithstanding  its  anomalous  RMS.  Model 
EOT1  la  may  be  the  slightly  preferred  model  for  these  56  test  sites,  since  both  its  median  differences  and  its 
RMS  differences  are  generally  very  good  across  most  constituents.  Flowever,  FES1 2  has  the  most  bold  entries 
for  the  RMS  differences,  and  TPX08  appears  best  for  median  differences  for  all  constituents  except  for  M2. 

5.  High-Latitude  Tides 

Regions  poleward  of  the  Topex/Poseidon  and  Jason  maximum  latitude  of  66°,  and  where  sea  ice  and  ice 
shelves  interfere  with  satellite  altimeter  detection  of  the  ocean  surface,  have  lower  densities  of  high-quality 
satellite  and  in  situ  data  for  validation  and  assimilation.  These  regions  also,  typically,  have  poor  bathymetric 
grids,  especially  for  regions  under  the  large  Antarctic  ice  shelves.  In  this  section  we  report  separately  on  tide 
model  accuracy  for  the  Arctic  and  Antarctic  regions  and  briefly  describe  future  prospects  for  improvement. 

5.1.  Arctic  Ocean 

Tide  models  were  tested  in  the  Arctic  Ocean  (latitudes  >65°N)  using  a  data  set  of  harmonic  constants  for 
M2,  S2,  K,,  and  O,,  determined  previously  from  240  tide  gauges,  with  the  number  of  stations  being  differ¬ 
ent  for  the  different  constituents.  (The  data  set,  developed  by  Z.  Kowalik  and  A.  Proshutinsky,  is  available 
at  http://www.ims.uaf.edu/tide).  Due  to  harsh  conditions,  Arctic  tide  gauges  are  often  placed  in  sheltered, 


M2  -  GOT4.8 

mean  =  4.25  cm 
median  =  2.45  cm 


T- 


X 


4- 


10  15  20  25  30  35  40 

Vector  difference  (cm) 

Figure  3.  Flistograms  of  M2  vector  differences  (absolute  magnitudes) 
between  195  shelf  tide  stations  and  one  of  the  more  accurate  shelf  tide 
models  (DTU10)  and  one  of  the  less  accurate  (GOT4.8)  according  to  the 
statistics  of  Table  4.  Both  models  evidently  are  reasonably  accurate  for 
most  stations,  with  median  vector  differences  less  than  2.5  cm.  Flowever, 
both  also  show  large  model-gauge  differences  for  a  handful  of  stations. 
Thus,  most  tested  models  are  reasonably  accurate  on  the  shelves  and 
in  shallow  seas,  but  large  errors  occur  in  a  few  locations,  which  are  not 
always  predictable  and  which  may  surprise  unwary  users. 
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Table  5.  Tide  Signal  and  RMS  Model  Differences  (cm)  With  56  Coastal  Tide  Gauges3 


Qi 

o, 

Pi 

K, 

n2 

m2 

s2 

K2 

RSS 

m4 

Signal 

2.0 

10.5 

5.3  17.1  12.6  60.5 

Root-Mean-Square  Differences 

11.2 

6.2 

1.9 

CSR3.0 

0.70 

2.77 

1.46 

5.03 

4.62 

20.74 

7.84 

2.31 

23.53 

GOT4.8 

0.46 

1.01 

0.57 

1.80 

1.92 

7.00 

3.53 

1.34 

8.46 

1.23 

OSU12 

1.14 

1.68 

1.20 

3.20 

2.15 

9.01 

5.31 

1.41 

11.48 

0.99 

DTU10 

0.62 

1.29 

0.73 

2.08 

1.72 

5.24 

2.68 

1.40 

6.82 

1.23 

EOTIIa 

0.54 

1.32 

0.85 

2.38 

1.78 

4.50 

2.84 

1.49 

6.49 

1.04 

HAM  12 

0.29 

1.42 

0.71 

2.65 

2.27 

14.63 

5.41 

1.80 

16.16 

FES12 

0.32 

0.89 

0.61 

1.65 

1.74 

6.60 

2.27 

0.77 

7.50 

1.49 

TPX08 

0.43 

1.13 

0.93 

2.01 

3.34 

15.65 

7.79 

2.12 

18.10 

1.68 

Bootstrap  a 

0.10 

0.23 

0.10  0.36  0.29  0.95 

Median  Absolute  Differences 

0.69 

0.32 

0.18 

CSR3.0 

0.25 

0.56 

0.32 

0.94 

0.63 

2.48 

1.31 

0.57 

3.15 

GOT4.8 

0.14 

0.42 

0.27 

0.62 

0.60 

1.30 

0.80 

0.37 

1.86 

0.21 

OSU12 

0.35 

0.46 

0.26 

0.62 

0.59 

1.19 

1.58 

0.34 

2.27 

0.23 

DTU10 

0.21 

0.37 

0.37 

0.67 

0.42 

1.28 

0.77 

0.42 

1.83 

0.20 

EOTIIa 

0.18 

0.36 

0.30 

0.61 

0.41 

1.06 

0.65 

0.32 

1.56 

0.25 

HAM  12 

0.12 

0.46 

0.27 

0.68 

0.40 

1.02 

0.66 

0.28 

1.57 

FES12 

0.18 

0.43 

0.44 

0.66 

0.47 

1.36 

0.66 

0.24 

1.85 

0.21 

TPX08 

0.13 

0.26 

0.25 

0.58 

0.28 

1.49 

0.49 

0.21 

1.75 

0.14 

Bootstrap  a 

0.02 

0.05 

0.03 

0.09 

0.06 

0.23 

0.11 

0.04 

0.03 

aModels  listed  in  order  of  increasing  spatial  resolution.  Bold  font  marks  models  within  Icr  of 
the  smallest  RMS,  for  a  particular  constituent. 


narrow  fjords,  thus  lowering  the  chance  that  these  locations  lie  within  model  domains.  For  the  investiga¬ 
tion  it  was  required  that  at  least  one  valid  model  point  was  available  close  to  the  tide  gauge  position  in 
order  to  interpolate  the  model  results  and  to  keep  the  tide  gauge  data,  and  because  of  this  step  many  tide 
gauges  had  to  be  dropped  from  further  analyses.  Because  the  original  number  of  tide  gauges  for  which 

individual  constituents  are  available  dif- 


Distance  from  tide  gauge  (km) 

Figure  4.  Amplitude  of  the  M2  tide  near  Duck  Pier,  North  Carolina, 
which  sits  on  a  barrier  island  with  the  Atlantic  Ocean  on  the  east 
(right)  and  Albemarle  Sound  on  the  west  (left).  The  tide  gauge 
amplitude  is  47.6  cm  and  marked  by  the  large  open  circle.  The 
gauge  accurately  measures  the  Atlantic  tide;  the  tide  in  the  sound 
is  much  smaller.  Colored  lines  display  model  amplitudes  east  and 
west  of  the  tide  gauge,  with  colored  dots  at  the  resolution  of  each 
model  (e.g.,  GOT4.8  plotted  every  0.5°);  unconnected  dots  imply 
that  one  or  more  model  values  are  missing  (i.e.,  are  categorized 
as  land).  Most  models  are  confused  by  the  presence  of  the  barrier 
island  and  disagree  sharply  in  Albemarle  Sound.  Models  DTU10  and 
EOT1 1  a  closely  follow  their  prior,  FES2004. 


fered,  the  final  number  of  stations  used 
in  the  comparison  between  constituents 
differs  as  well. 

Another  source  of  error  is  a  seasonal  varia¬ 
tion  of  polar  tides  inferred  by  the  seasonally 
varying  ice  coverage.  In  Arctic  regions 
this  effect  causes  changes  in  tidal  ampli¬ 
tudes  of  O  (1-1 0  cm).  Since  the  modern 
tide  models  represent  an  annual  mean 
state  but  most  of  the  observation  do  not 
account  for  the  observed  seasonality,  this 
implies  some  inconsistency  in  this  model- 
observation  evaluation. 

Because  the  Arctic  tide  gauge  data  set  has 
not  undergone  extensive  quality  control, 
we  edited  the  data  prior  to  the  model  eval¬ 
uation.  This  was  performed  on  the  reduced 
60  tide  gauge  set  common  to  the  five  mod¬ 
els  HAM12,  FES12,  EOT1  la,  DTU10,  and 
TPX08;  in  cases  where  more  than  three 
models  showed  RSS  differences  for  the  four 
largest  constituents  greater  than  1 5  cm  with 
respect  to  the  gauge  data,  the  station  was 
removed  from  the  analysis. 
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Figure  5.  Magnitudes  of  vector  differences  (in  cm)  between  in  situ  and  bilinearly  interpolated  M2  tide  amplitude  for  60  tide  gauges  for  (a)  HAM12,  (b)  FES12, 
(c)  EOT1  la,  (d)  DTU10,  and  (e)  TPX08  over  the  Arctic.  Only  20  gauges  could  be  used  for  the  interpolation  for  (f)  GOT4.8.  OSU12  is  identical  to  GOT4.8  and  is 
therefore  not  shown. 


This  test  eliminated  the  gauges  Kanin  and  Bugrino  (both  Russia)  and  Iginiarfic  (West  Greenland).  Figure  5 
shows  the  resulting  final  vector  differences  between  in  situ  and  predicted  M2  elevation  amplitude  for  the 
common  tide  gauges.  Due  to  a  restrictive  model  domain  for  the  GOT4.8  model,  only  20  gauges  were  used 
for  the  comparison  for  this  model;  and,  as  OSU1 2  defaults  to  GOT4.8  at  high  latitude,  only  the  results  for 
GOT4.8  are  shown  (Figure  5f).  Figure  5e  illustrates  that  TPX08  generally  has  very  small  differences  with  the 
tide  gauges,  a  consequence  of  several  of  these  being  assimilated  into  the  model.  Flowever,  a  few  gauge 
positions  show  enhanced  differences  in  the  southern  part  of  the  Barents  Sea,  e.g.,  Vardo  (Norway),  which 
has  a  very  large  difference  for  this  model  only,  pointing  toward  a  problem  with  the  gauges  being  assimilated 
into  the  model  in  this  region.  In  addition,  FIAM1 2  and  FES12  both  show  poor  agreement  with  three  gauges 
in  Nares  Strait  at  the  northern  end  of  Baffin  Bay.  (Locations  are  shown  by  numbers  1, 3,  and  4  in  Figure  5b.) 

The  RMS  and  RSS  for  the  four  main  tide  constituents  over  the  Arctic  are  shown  in  Table  6  for  the  reduced 
set  of  gauges  common  to  all  seven  models.  As  can  be  seen  from  the  table,  DTU1 0  shows  the  best  overall 
agreement  with  tide  gauge  data  with  an  RSS  of  5.7  cm,  whereas  all  other  models  have  RSS  values  in  the 
range  of  6-8  cm.  It  should  also  be  noted  that  all  gauges  are  coastal  gauges  and  consequently  the  numbers 
should  be  compared  with  previous  comparisons  with  coastal  gauges  for  lower  latitudes  (see  section  4.3). 

For  the  four  models  not  assimilating  the  tide  gauge  information  (HAM12,  FES12,  EOT1  la,  and  DTU10),  the 
mean  RMS  difference  amounts  to  16%  of  the  M2  average  amplitude  and  more  than  30%  for  S2,  K,,  and  0,. 
This  illustrates  the  problem  that  altimetric  tidal  models  face  in  the  Arctic.  Only  Sun-synchronous  satellites 
(ERS-1,  ERS-2,  and  Envisat),  for  which  the  S2  is  phase  locked  and  aliases  into  a  constant,  have  been  available 
in  the  Arctic.  Similarly,  K,  and  P,  are  aliased  into  the  annual  signal. 

Tide  models  were  also  tested  against  three  bottom  pressure  time  series  in  the  Arctic  close  to  the  geographic 
North  Pole:  ABPR1  at  89.2543°N,  60.3597°E;  ABPR3  at89.2475°N,  148.1257°E;  and  ABPR4  at  89.9767°N, 
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Table  6.  RMS  and  RSS  Differences  (Both  in  cm)  Between 
the  Interpolated  Tidal  Signal  and  Common  Tide  Gauges 
for  the  Four  Major  Arctic  Ocean  Tidal  Constituents3 


No.  of  Stations 


0, 

13 


K, 

16 


M2 

20 


S2 

15 


Signal 

4.81 

10.35 

30.46 

11.45 

RSS 

GOT4.8 

1.78 

3.20 

4.83 

2.24 

6.47 

OSU12 

1.77 

3.17 

4.82 

2.22 

6.43 

DTU10 

1.60 

2.89 

3.91 

2.56 

5.72 

EOTIIa 

1.36 

2.85 

4.61 

2.44 

6.09 

HAM12 

1.25 

3.00 

4.11 

7.11 

8.83 

FES12 

1.52 

3.30 

4.66 

3.36 

6.80 

TPX08b 

1.19 

1.44 

5.89 

1.93 

6.47 

aThe  average  amplitude  (in  cm)  for  each  tidal  con¬ 
stituents  of  the  Arctic  is  listed.  The  values  are  given 
for  the  largest  common  set  of  tide  gauges  to  five  of 
the  models. 

bTPX08  model  assimilates  a  subset  of  the  selected  in 
situ  measurements. 


1 78.23° E.  For  each  bottom  pressure  gauge,  more 
than  1  year  of  data  was  used  to  derive  the  tidal 
constituents.  Since  the  tidal  signal  is  relatively 
small  at  the  North  Pole  (M2  amplitude  varies 
between  5.3  and  6.1  cm  among  the  gauges),  the 
RMS  differences  and  RSS  differences  between 
the  interpolated  tidal  signal  and  bottom  pres¬ 
sure  gauges  for  the  eight  major  tidal  constituents 
range  only  from  1 .2  to  2.3  cm,  with  HAM1 2  and 
FES1 2  having  the  smallest  RSS.  All  other  models 
have  higher  RSS  due  to  large  differences  for  indi¬ 
vidual  constituents  (K,  for  GOT4.8,  OSU12,  and 
TPX08,  and  S2forDTU10). 

5.2.  Antarctic  Seas 

A  total  of  1 02  in  situ  measurements  are  available 
from  seas  surrounding  Antarctica.  Data  sources 
include  bottom  pressure  recorders;  coastal  tide 
gauges;  and  GPS,  gravimeter,  and  tiltmeter  mea¬ 


surements  on  the  floating  ice  shelves  [King  and Padman,  2005;  King  etai.,  201  la,  201 1  b].  All  derived  tidal 
constituents,  along  with  information  about  record  length,  measure  type,  and  references  are  included  in 
the  Antarctic  Tide  Gauge  Database  (http://www.esr.org/antarctic_tg_index.html).  The  gravimeter  and  tilt- 
meter  records  are  expected  to  be  less  reliable,  generally  being  based  on  shorter  and  older  records  and 
a  greater  number  of  uncertain  corrections  used  in  the  reduction.  Unfortunately,  the  gravimeter  records 
make  up  a  large  fraction  of  the  Ross  Ice  Shelf  records.  Flowever,  the  tide  gauges  are  generally  quality  con¬ 
trolled.  To  verify  this,  the  identical  comparison  between  models  and  gauges  to  that  for  the  Arctic  Ocean 
was  performed. 


The  maps  of  vector  differences  between  interpolated  and  observed  M2  signals  (Figure  6)  suggest  that,  espe¬ 
cially  in  regions  covered  by  the  large  ice  shelves,  all  models  have  problems  in  producing  the  correct  tidal 
signal.  All  the  models  show  vector  differences  larger  than  1 0  cm  on  the  Filchner-Ronne  Ice  Shelf,  where  the 
amplitude  of  M2  exceeds  1  m.  Similarly,  on  the  Ross  Ice  Shelf,  all  the  models  except  FIAM1 2  and  TPX08  have 
vector  differences  for  M2  exceeding  1 0  cm,  representing  a  substantial  relative  error  as  the  amplitude  is  less 
than  30  cm.  Flowever,  the  gravimetry  data  here  are  less  reliable. 


The  calculated  RMS  differences  between  in  situ  measurements  and  model  predictions  are  listed  in  Table  7. 
Because  OSU1 2  is  not  available  on  the  ice  shelves,  it  is  excluded  from  this  comparison.  The  RMS  differ¬ 
ences  for  the  full  set  of  common  gauges  are  3-4  cm,  2-5  cm,  2-6  cm,  and  1-5  cm  for  the  M2,  S2,  K,,  and  01 
constituents,  respectively.  Most  models  have  large  RSS  on  the  Larsen  C  and  Filchner-Ronne  ice  shelves  in 
the  Weddell  Sea.  TPX08  exhibits  the  lowest  RMS  for  most  constituents  and  clearly  the  lowest  overall  RSS, 
consistent  with  its  assimilation  of  many  of  these  in  situ  data  records. 


Table  7  also  shows  the  average  constituent  amplitude  over  the  tide  gauges.  In  contrast  to  the  Arctic  Ocean, 
amplitudes  of  the  K,  and  O,  constituents  around  sections  of  Antarctica  can  be  large;  the  four  major  con¬ 
stituents,  M2,  S2,  K,,  and  01  display  roughly  similar  amplitudes  around  Antarctica  (although  this  assessment 
is  influenced  by  the  irregular  distribution  of  the  in  situ  data).  For  the  five  models  that  do  not  assimilate  data, 
i.e.,  FI  AM  12,  FES12,  EOT1  la,  DTU10,  and  GOT4.8,  the  mean  RMS  error  for  M2  is  17%  of  the  M2  average  ampli¬ 
tude.  This  RMS  error  increases  to  more  than  20%  for  most  other  constituents,  with  a  value  of  nearly  30%  for 
the  smaller  N2  tide. 

To  illustrate  the  magnitude  of  the  misfit,  Figure  7  compares  GPS-based  elevations  [King  et  at.,  2011a,  201 1  b] 
with  model  time  series  for  the  two  sites  FR03  and  FR09  located  on  the  Filchner-Ronne  Ice  Shelf  in  the 
Weddell  Sea.  Ocean  tide  loading  displacement  and  atmospheric  pressure  ("inverse  barometer")  effects  have 
been  removed  from  the  GPS  time  series,  leaving  a  relatively  smooth  time  series  dominated  by  ocean  tides. 
Tidal-band  differences  of  several  tens  of  centimeters  are  evident  over  the  time  window  and  vary  substan¬ 
tially  in  amplitude  and  phase  between  models;  such  differences,  and  larger  ones,  are  common  to  all  GPS 
sites  on  the  Larsen  C  and  Filchner-Ronne  ice  shelves.  In  contrast,  Kim  et  al.  [201 1]  and  Oreiro  et  al.  [2014] 
report  model-data  differences  for  other  Antarctic  regions  that  are  much  smaller  than  those  shown  here, 
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Figure  6.  Magnitudes  of  vector  differences  (in  cm)  between  in  situ  and  interpolated  M2  tide  amplitude  for  86  sites  that  could  be  used  by  all  models  except  for  the 
OSU12,  which  excludes  Antarctic  ice  shelves  from  its  domain  and  therefore  is  not  shown. 

although  some  of  the  models  they  test  are  the  same  as  those  that  we  use  here.  These  results  highlight  large 
spatial  differences  in  tide  model  accuracy,  in  part  depending  on  distance  from  altimeter-based  model  con¬ 
straints,  complexity  of  the  local  coastline,  and  water  depth  errors  especially  under  large  ice  shelves.  We  note 
that  King  et  al.  [201 1  a]  investigated  seasonal  variation  in  the  amplitudes  of  M2  and  O,  in  the  Weddell  Sea 
and  found  them  to  be  negligible. 

6.  Tests  of  Modern  Tide  Models  Against  Satellite  Data 

Tide  model  tests  against  in  situ  data  are  primarily  local  tests  that  could  be  hampered  by  details  of  the  bot¬ 
tom  topography  in  the  models  or  other  uncertain  local  model  parameters.  Such  tests  can  also  be  limited  by 
the  quality  or  sparseness  of  the  in  situ  data,  a  particular  concern  in  polar  regions.  Testing  the  impact  of  tide 
models  against  various  kinds  of  satellite  measurements  provides  some  crucial  complementary  assessments. 
Some  of  these  tests  (e.g.,  the  effect  on  satellite  orbit  perturbations)  are  limited  to  an  indirect,  integral,  quality 
assessment,  but  these  can  still  be  useful  and  are  in  any  event  of  interest  to  the  space  geodetic  community. 
Such  tests  are  reported  in  this  section. 

6.1.  Satellite  Tracking  Tests 

The  orbits  of  artificial  satellites  are  significantly  perturbed  by  the  gravitational  effects  of  ocean  tides  [e.g., 
Lambeck,  1988].  Monitoring  such  orbital  perturbations  with  precise  satellite  tracking  measurements  thus 
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Table  7.  RMS  and  RSS  Differences  (Both  in  cm)  Between  the  Interpolated  Antarctic  Tidal 
Signal  and  Common  Tide  Gauges3 


No.  of  Stations 
Signal 

Qi 

36 

6.79 

o, 

49 

30.50 

Pi 

37 

10.80 

Ki 

48 

32.43 

n2 

37 

8.12 

m2 

49 

37.62 

s2 

48 

27.65 

k2 

36 

9.28 

RSS 

GOT4.8 

0.88 

3.52 

1.04 

4.01 

0.94 

3.80 

3.49 

0.74 

7.64 

DTU10 

1.09 

4.70 

1.87 

5.38 

1.70 

3.10 

4.74 

1.35 

9.62 

EOTIIa 

1.23 

4.92 

2.05 

5.83 

1.82 

3.48 

4.80 

1.36 

10.21 

HAM12 

0.78 

2.55 

1.37 

3.36 

1.06 

3.20 

3.26 

0.96 

6.58 

FES12 

1.50 

4.55 

1.78 

4.63 

1.06 

4.27 

4.53 

1.13 

9.42 

TPX08b 

0.53 

1.48 

0.74 

2.22 

0.76 

3.14 

2.37 

0.49 

4.92 

3The  average  amplitude  (in  cm)  for  each  tidal  constituents  of  the  Antarctic  is  listed. 
The  values  are  given  for  the  largest  common  set  of  tide  gauges  to  five  of  the  models 
(excluding  TPX08). 

bDenotes  a  model  that  assimilates  a  subset  of  the  selected  in  situ  measurements. 


provides  an  independent  method  to  assess  accuracies  of  global  tidal  models.  The  following  subsections 
describe  two  complementary  types  of  tracking:  ground  to  satellite  and  satellite  to  satellite.  For  the  former 
we  use  satellite  laser  ranging  (SLR)  from  ground  observatories;  for  the  latter  we  use  microwave  range  rates 
between  the  two  Gravity  Recovery  and  Climate  Experiment  (GRACE)  satellites.  Owing  to  the  inherent  limita¬ 
tions  of  spaceborne  gravity  measurements,  both  are  sensitive  only  to  long-wavelength  components  of  tidal 
mass  variations.  The  GRACE  satellites  are  sensitive  primarily  to  scales  above  400  km  [Rowlands  et  ai,  2005]. 
SLR  is  sensitive  to  only  a  few  low-degree  spherical  harmonic  components,  generally  of  order  1  for  diurnal 
constituents  and  order  2  for  semidiurnal  constituents. 

To  analyze  the  tidal  residuals  in  either  type  of  satellite  tracking  data  requires  gravitational  models  of  nontidal 
variability  which  are  as  complete  as  possible.  For  the  tests  below,  similar  background  models  were  used, 
including  GRACE-based  models  of  the  static  and  time-variable  gravity  field  and  models  of  rapid  atmospheric 
variability  implied  by  operational  models  of  the  ECMWF.  The  GRACE  tests  also  used  models  of  ocean  and 
continental  hydrological  variability,  and  they  used  higher-degree  expansions  for  all  fields.  In  all  cases  the 
body  tides  were  modeled  following  2010  Conventions  of  the  International  Earth  Rotation  and  Reference 
Systems  Service. 

Each  tested  tide  model  was  expanded  into  spherical  harmonic  series.  Only  the  five  major  constituents 
M2,  S2,  N2,  O,,  and  K,  were  used  for  each  model,  with  25  other  minor  tidal  constituents  (or  in  some  cases, 
tidal  groups)  based  on  GOT4.8,  except  for  long-period  tides  which  were  based  on  either  TPX06.2  or  on 


FR03  IBE  and  OTL  removed  from  GPS 


FR09  IBE  and  OTL  removed  from  GPS 


-0.4  - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

40  41  42  43  44  45  46  47  48 

Day  of  Year  (2008) 


Figure  7.  Comparison  of  GPS-based  elevations  [ King  et  al.,  2011a,  201 1  b]  with  model  time  series  for  the  two  sites  FR03  and  FR09  located  on  the  Filchner-Ronne 
Ice  Shelf  in  the  Weddell  Sea. 
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self-consistent  equilibrium  models.  Nodal  and  perigee  modulations  of  lunar  tides  were  accounted  for  by 
standard  methods.  For  GRACE  the  degree  of  expansion  was  taken  to  90;  for  the  SLR  satellites,  to  20.  The 
computation  of  the  gravitational  potential  from  the  ocean  tide  expansions  accounts  for  crustal  loading 
[e.g.,  Lambeck,  1 988,  equation  (6.1 .24)],  requiring  an  adopted  set  of  loading  Love  numbers  k' .  The  same 
computation  requires  knowledge  of  seawater  density,  which  we  took  as  a  constant  of  p  =  1031  kg  m~3. 

6.1.1.  Satellite  Laser  Ranging  Tests 

The  impact  of  individual  tide  models  on  satellite  orbit  computations  was  tested  against  precise  SLR  tracking 
data  from  four  geodetic  satellites  orbiting  at  a  variety  of  altitudes  and  inclinations:  Lageos,  Starlette,  Stella, 
and  Larets.  All  four  satellites  are  heavy,  "cannonball,"  geodetic  satellites  covered  with  cube  corner  reflectors. 
For  all  four  satellites  we  undertook  a  consistent  set  of  analyses  based  on  261  seven-day  orbital  arcs  over  the 
period  2004-2008.  Many  SLR  analyses  estimate  empirical  acceleration  adjustments  at  resonant  frequencies 
(e.g.,  once  per  revolution)  to  accommodate  force  model  errors;  however,  these  kinds  of  adjustments  poten¬ 
tially  dilute  the  sensitivity  to  tidal  perturbations  in  the  orbit,  so  we  here  avoided  such  adjustments.  The 
metric  of  this  test  is  to  identify  the  tide  model  that  produces  orbits  with  the  smallest  SLR  tracking  residuals. 

Because  SLR  tracking  is  not  continuous  in  time  but  rather  relies  on  occasional  satellite  overpasses  of  the 
laser  observatories,  the  SLR  data  are  sensitive  essentially  to  only  long-period  perturbations  in  satellite  orbits. 
These  perturbations  depend  on  the  orbital  elements  of  a  satellite,  so  each  satellite  is  more  sensitive  to 
some  tidal  constituents  than  others.  We  therefore  should  not  expect  results  to  be  completely  consistent 
among  the  four  satellite  tests.  For  example,  Starlette,  specially  designed  for  sensing  tidal  perturbations,  is 
most  sensitive  to  K,,  and  then  (in  order)  S2,  M2,  and  P,,  with  the  K,  perturbations  having  periods  near 
90  days  [Williamson  and  Marsh,  1985].  Lageos  is  also  sensitive  to  K1(  K2,  and  S2,  but  the  K,  periodicities  are 
near  3  years  and  therefore  potentially  confounded  by  nontidal  mass  variations  [Lambeck,  1988;  Tables  6.8]. 
Lageos  is  also  at  far  higher  altitude  than  the  others  and  is  therefore  sensitive  to  only  the  very  lowest  degree 
gravitational  terms.  Both  Stella  and  Larets  are  in  near-Sun-synchronous  orbits. 

Table  8  presents  the  mean  RMS  residuals  in  SLR  tracking  data  for  seven  tested  tide  models,  all  of  them 
modern  altimetric-based  models  except — to  assess  historical  progress — the  old  NSWC  [Schwiderski,  1979] 
model.  Each  entry  in  the  table  gives  the  mean  RMS  over  the  261  arcs,  with  the  quoted  uncertainties 
being  the  standard  errors  in  those  means.  For  each  satellite  the  NSWC  model  clearly,  and  reassuringly, 
yields  the  largest  SLR  residuals.  As  an  overall  assessment  of  the  SLR  tests,  Table  8  suggests  that  models 
EOT1  la,  FIAM12,  and  GOT4.8  all  perform  well  for  all  four  satellites.  TPX08  is  especially  good  for  Starlette  but 
appears  slightly  less  accurate  for  the  two  Sun-synchronous  satellites.  The  OSU12  model  yields  less  accurate 
computed  orbits,  probably  because  it  lacks  tidal  data  under  Antarctic  ice  shelves  (cf.  results  in  next  section). 

6.1.2.  GRACE  Tests 

Satellite-to-satellite  ranging  data  in  GRACE  are  continuous  in  time,  so  (unlike  the  SLR  test)  GRACE  data  are 
sensitive  to  both  short-period  and  long-period  perturbations  in  the  orbits.  For  a  given  mass  anomaly  of 
sufficient  size,  GRACE  perturbations  can  be  localized  over  that  mass,  which  of  course  is  the  source  of  its  great 
value  for  monitoring  mass  changes  [Tapley  et  al.,  2004].  Our  GRACE  tests  are  based  on  the  semiqualitative 
approach  described  by  Ray  et  al.  [2009].  Using  similar,  but  more  complete,  background  models  as  described 
above  for  the  SLR  tests,  with  higher-degree  tide  expansions,  we  have  computed  GRACE  range  rate  residuals 
for  the  7year  period  2004-2010.  The  processing  of  the  range  rate  data  follows  steps  described  by  Luthcke  et 
al.  [2006], 

The  entire  processing  was  repeated  for  each  tested  tide  model.  These  range  rate  residuals  were  then  differ¬ 
entiated  along  track  and  low-pass  filtered  to  generate  range  acceleration  residuals,  and  these  residuals  were 
then  binned  into  small  geographical  areas  and  tidally  analyzed — see  Ray  et  al.  [2009]  for  further  details. 
Regions  displaying  significant  tidal  energy  in  the  residuals  suggest  regions  of  likely  errors  in  the  tested 
model.  The  test,  however,  is  qualitative  because  the  final  maps  are  still  in  terms  of  GRACE  measurement 
residuals  and  not  in  terms  of  tidal  elevations  of  the  sea  surface,  which  would  require  an  additional  inversion 
to  determine  (for  recent  attempts  at  such  inversions,  see  Han  et  al.  [2007]  and  Killett  et  al.  [201 1]).  Moreover, 
the  signature  of  a  mass  anomaly  in  range  acceleration  includes  significant  sidelobes  [e.g.,  Ray  et  al.,  2009, 
Figure  1],  which  act  to  extend  anomalous  regions  over  wider  areas,  including  over  land;  this  point  must  be 
kept  in  mind  when  interpreting  our  maps. 

Tidal  amplitudes  of  GRACE  range  acceleration  in  units  of  nm  s-2  are  shown  in  Figures  8-1 1  for  constituents 
M2,  K1(  O,,  and  S2,  and  the  global  means  of  these  residual  amplitudes  are  tabulated  in  Table  9.  Figure  8  (M2) 
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Table  8.  Mean  RMS  (cm)  of  Satellite  Laser  Ranging  Residuals  to  Four 
Geodetic  Satellites3 


Lageos-1 

Stella 

Starlette 

La  rets 

NSWC 

1.447  ±0.021 

4.81 9  ±0.011 

5.1 79  ±0.009 

5.040  ±0.011 

GOT4.8 

1.231  ±0.020 

2.875  ±  0.007 

3.043  ±  0.006 

3.263  ±0.009 

OSU12 

1.259  ±0.020 

3.714±0.009 

3.675  ±0.008 

3.798  ±0.010 

EOTIIa 

1.235  ±0.020 

2.991  ±0.007 

3.1 35  ±0.007 

3.341  ±0.009 

HAM  12 

1.229  ±0.020 

3.025  ±0.007 

3.061  ±0.006 

3.268  ±0.009 

FES12 

1.238  ±0.020 

3.203  ±0.008 

3.1 77  ±0.007 

3.449  ±0.010 

TPX08 

1.227  ±0.020 

3.1 57  ±0.008 

2.993  ±  0.006 

3.599  ±0.011 

tabulated  RMS  for  each  satellite  is  the  mean  of  the  RMS  residuals  com¬ 
puted  for  each  of  261  seven-day  arcs  over  the  period  2004-2008.  Quoted 
uncertainties  are  standard  errors  in  these  means.  Bold  font  marks  models 
within  1  standard  error  of  the  smallest  RMS,  for  each  constituent.  DTU10 
model  was  not  used  in  this  comparison. 


shows  the  expected  high-latitude  problems,  where  the  tide  models  are  known  to  be  less  accurate  owing  to 
the  lack  of  satellite  altimetry.  In  addition,  however,  the  resolution  of  GRACE  is  sufficiently  high  to  allow  us  to 
resolve  relatively  isolated,  near-coastal  anomalies  in  lower  latitudes.  For  example,  EOT1  la  displays  small  M2 
anomalies  along  western  India,  whereas  HAM12  (which  assimilated  the  EOT1  la  tides)  has  mostly  corrected 
these  problems.  Most  models  (except  TPX08)  show  small  anomalies  along  the  coast  of  British  Columbia. 
These  figures  can  therefore  guide  modelers  to  locations  where  improvements  may  be  needed. 

Nevertheless,  it  is  also  possible  that  some  anomalies  shown  in  Figure  8  reflect  GRACE  (or  GRACE  processing) 
errors  and  not  tide  model  errors.  For  example,  all  six  models  show  anomalies  near  New  Zealand,  where  the 
M2  wave  is  large  and  of  relatively  short  wavelength  as  it  rotates  around  the  islands;  our  use  of  a  spherical 
harmonic  series  limited  to  degree  90  may  be  inadequate  to  capture  the  full  M2  signal  near  New  Zealand.  The 
large  anomalous  regions  of  the  North  Atlantic  may  also — given  that  all  six  models  are  similar — represent 
systematic  errors  in  GRACE  processing.  For  example,  errors  in  the  adopted  density  of  seawater  in  regions  of 
large-amplitude  tides,  such  as  off  western  Europe,  may  be  partly  to  blame — see  Figure  10  of  Rayetal.  [2009]. 

Figure  9  emphasizes  that  all  models  have  significant  errors  around  Antarctica.  The  extent  of  the  errors  in 
TPX08  may  be  somewhat  misleading,  owing  to  sidelobes  of  the  main  anomaly,  especially  in  the  Ross  Sea. 
Sidelobes  are  also  suggested  in  the  FIAM1 2  panel  in  the  western  Ross  Sea.  Because  EOT1  la  defaults  to 
FES2004  in  high  latitudes,  Figure  9  suggests  that  FES2004  was  more  accurate  than  FES1 2  in  the  Arctic  K, . 
Similarly,  Figure  10  (O,)  suggests  that  FES2004  is  more  accurate  than  FES1 2  in  the  Weddell  Sea,  but  the 
reverse  is  true  in  the  Ross  Sea,  where  FES12  appears  to  be  the  most  accurate  of  all  six  models. 


EOT  1 1 


GOT4 


HAM  12 


OSU12 


M2  !  !  !  !  !  !  ! 
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Figure  8.  Amplitudes  at  the  M2  tidal  frequency  in  7  years  of  GRACE  range  acceleration  residuals,  based  on  six  different 
tidal  models.  Locations  showing  significant  amplitudes  suggest  errors  in  the  corresponding  tide  model. 
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The  OSU12  model  performs  noticeably  poorly  around 
Antarctica  in  all  constituents  for  one  obvious  reason:  the 
model  does  not  include  regions  under  the  permanent 
Antarctic  ice  shelves.  Tides  in  those  regions  are  generally 
very  large  (although  M2  is  small  under  the  Ross  Ice  Shelf), 
so  neglecting  these  regions  will  induce  significant  residuals 
in  GRACE  ranging  data.  The  same  error  has  also  undoubt¬ 
edly  contributed  to  OSU1 2's  somewhat  higher  SLR  residuals 
shown  previously  in  Table  8. 

In  S2,  shown  in  Figure  11,  the  overall  background  of  the 
residuals  is  high,  especially  when  compared  with  the  very 
quiet  O,  residuals.  The  noise  is  likely  induced  by  various  systematic  errors  related  to  the  GRACE  orbit  plane 
alignment  with  the  Sun  (so-called  /?'  errors);  the  orbit  plane  precession  period  is  identical  to  the  S2  alias 
period.  In  addition,  in  low  latitudes,  errors  in  the  S2  atmospheric  tide  contribute  but  only  to  the  extent  that 
errors  are  identical  in  all  six  panels.  In  fact,  altimeter-based  models  themselves  tend  to  display  /?'-like  errors; 
Ray  [2013]  discusses  several  of  these,  including  errors  in  most  models  from  inadequate  air  tide  modeling  in 
theT/P  dry-troposphere  correction  and  the  lack  in  all  models  of  a  correction  for  crustal  loading  by  the  S2  air 
tide.  (We  note  that  a  small  anomaly  in  all  panels  over  the  Amazon  basin  is  caused  by  an  error  in  the  adopted 
hydrological  model  which  occurs  at  the  161  day  S2  alias  period.) 

Both  Table  9  and  Figures  8-1 1  suggest  that  the  FIAM12  model  may  be  the  best  for  processing  GRACE  range 
rate  data,  with  EOT1 1  a  and  GOT4.8  close  behind.  We  note  that  TPX08,  which  appears  perhaps  the  best  in 
lower  latitudes,  has  relatively  large  errors  over  the  Weddell  and  Ross  Seas;  this  is  surprising  since  TPX08  is 
the  only  model  that  assimilated  a  large  number  of  in  situ  data  in  those  regions.  It  is  possible  that  some  of 
those  in  situ  data,  such  as  old  gravimeter  stations  on  floating  ice,  may  be  problematic  because  of  errors  in 
ocean  loading  [King  and  Padman,  2005]  or  errors  in  instrument  tilt.  The  new  GPS  data  reported  by  King  et 
al.  [201  la,  201 1  b],  including  those  data  shown  in  Figure  7,  were  not  assimilated  into  TPX08  (or  any  other 
model),  so  future  improvements  in  some  of  these  problematic  regions  highlighted  by  GRACE  should  be 
forthcoming. 

6.2.  Satellite  Altimeter  Tests 

Variance  reduction  tests  on  altimetric  sea  level  anomalies  constitute  a  standard  method  for  evaluating 
ocean  tide  models.  Tests  of  our  seven  primary  global  models  are  discussed  in  this  section. 

Each  tested  ocean  tide  model  was  separately  employed  to  generate  corrections  to  sea  level  anomalies 
from  a  number  of  different  satellite  missions.  In  each  case  we  have  used  the  eight  major  tidal  constituents 
(M2,  S2,  N2,  K2,  K, ,  O, ,  P-, ,  and  Q, ),  with  all  long-period  constituents  held  fixed  and  other  constituents  (e.g., 


Table  9.  Mean  GRACE  Range  Acceler¬ 
ation  Residual  Amplitudes  (nms-2)  at 
Tidal  Frequencies 


m2 

s2 

o. 

K, 

NSWC 

1.237 

0.610 

0.269 

0.426 

GOT4.8 

0.198 

0.161 

0.090 

0.135 

OSU12 

0.267 

0.190 

0.168 

0.197 

EOTIIa 

0.194 

0.161 

0.079 

0.132 

HAM12 

0.195 

0.165 

0.032 

0.146 

FES12 

0.206 

0.192 

0.092 

0.172 

TPX08 

0.212 

0.166 

0.107 

0.172 

GOT4 


HAM12 


OSU12 


K1 

0.00  0.25  0.50  0.75  1.00  1.25  1.50  nm/s2 

Figure  9.  As  in  Figure  8  but  for  the  K,  frequency. 
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Figure  10.  As  in  Figure  8  but  for  the  C^  frequency. 


TPX08 


MS4)  ignored.  This  is  admittedly  not  an  optimum  approach,  because  it  ignores  the  compound  tides  that  are 
included  with  some  of  our  tested  models  (e.g.,  FES201 2  is  distributed  with  seven  compound  tides).  Thus,  our 
tests  of  these  models  in  shallow  water  will  yield  residual  variances  somewhat  inflated  from  what  the  com¬ 
plete  models  would  presumably  have  yielded.  However,  the  approach  does  allow  consistent  testing  of  the 
major  constituents  with  identical  software. 

In  addition  to  the  tide  corrections  the  altimetric  data  were  processed  with  standard  sets  of  media,  instru¬ 
ment,  and  geophysical  corrections.  These  include  a  dynamic  atmospheric  loading  correction  via  the 
MoG2D  (2-D  Gravity  Wave  Model)  [Carrere  and  Lyard,  2003],  An  equilibrium  model  was  used  to  correct  for 
long-period  tides,  and  the  model  of  Matsumoto  et  at.  [2000]  was  used  for  all  load  tides.  Further  detailed 
comparison  methodologies  are  described  by  Fok  [2012]. 

The  satellite  altimeter  data  used  for  these  tests  include  Topex/Poseidon  (1 0  day  cycles  4-364),  Topex  Inter¬ 
leaved  Mission  (cycles  369-481),  Jason-1  (cycles  1-259),  Jason-1  Interleaved  Mission  (cycles  262-374), 
Jason-2  (cycles  1  -1 84),  and  the  Jason-1  Geodetic  Mission  (cycles  382-409);  GFO  (1 7  day  cycles  37-223), 
Envisat  (35  day  cycles  1 0-94),  and  the  Cryosat-2  (29  day  cycles  1  -40).  Much  of  these  data  were,  of  course, 
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Figure  11.  As  in  Figure  8  but  for  the  S2  frequency.  The  overall  inflation  of  background  noise  stems  from  a  multitude  of 
possible  //'-like  errors  in  both  GRACE  data  and  in  the  ocean  tide  models  themselves.  Some  signals,  however,  stem  from 
inadequate  modeling  of  hydrology,  such  as  the  small  anomalies  evident  in  the  Amazon  basin  which  may  exist  at  the 
160  day  S2  alias  period  for  GRACE. 
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Table  10.  Standard  Deviations  of  Residual  SSH  Anomaly  (in  cm)  Along  Satellite  Tracks  in  Shallow  Oceans  (Depth 
<  1000  m)  After  Applying  Ocean  Tide  Model  Corrections3 


Model 

T/P 

Interlaced 

GFO 

Envisat 

Jason-1 

Interlaced 

Jason-1 

Jason-2 

Jason-1 

GM 

CryoSat-2 

GOT4.8 

10.66 

10.37 

11.40 

10.55 

11.24 

10.99 

10.69 

10.67 

11.35 

DTU10 

10.66 

10.34 

11.45 

10.51 

11.18 

10.95 

10.67 

10.67 

11.30 

OSU12 

10.83 

10.12 

11.44 

10.65 

11.22 

10.98 

10.69 

10.72 

11.34 

EOTlIa 

10.82 

10.50 

11.65 

10.62 

11.37 

11.20 

10.78 

10.75 

11.44 

HAM12 

10.72 

10.47 

11.72 

10.76 

11.34 

11.15 

10.68 

10.67 

11.44 

FES12 

10.69 

10.34 

11.52 

10.59 

11.09 

10.96 

10.62 

10.61 

11.29 

TPX08 

10.82 

10.54 

11.74 

10.84 

11.44 

11.17 

10.78 

10.77 

11.58 

Before 

55.06 

54.82 

53.24 

49.49 

56.23 

56.38 

55.42 

54.63 

49.99 

aThe  bottom  line  is  the  respective  SSH  standard  deviation  before  ocean  tide  correction  was  applied. 


employed  in  the  development  of  the  models  but  not  all.  For  example,  all  models  are  independent  of  the 
Jason-1  geodetic  mission  data  and  all  are  independent  of  Cryosat-2  data.  All  the  altimeter  data  used  in  this 
section  were  extracted  from  the  Radar  Altimeter  Database  System  [Scharroo,  2008]. 

The  overall  statistics  for  the  altimetry  sea  level  anomaly  tests  as  a  function  of  ocean  tide  models  are  sum¬ 
marized  in  Tables  1 0  and  1 1  for  shallow  (<  1 000  m  water  depth)  and  deep-ocean  areas,  respectively.  One 
immediately  notices  that  differences  in  variance  among  the  different  missions  is  generally  greater  than  the 
differences  among  models  within  a  mission.  This  likely  arises  because  the  various  missions  sampled  the 
oceans  at  different  times,  and  it  also  may  reflect  differences  in  data  quality.  In  the  deep  ocean  the  resid¬ 
ual  altimeter  variances  are  nearly  indistinguishable  between  different  tide  models,  although  the  OSU12 
model  is  sometimes  anomalous:  it  displays  slightly  better  statistics  on  T/P-interlaced,  GFO,  and  Jason-1  mis¬ 
sions  but  slightly  weaker  statistics  for  other  missions,  while  other  models  show  similar  values  for  all  missions 
(within  0.05  cm). 

Considering  shallow-water  regions,  one  can  observe  some  differences  between  each  tidal  correction, 
although  differences  are  still  very  small:  FES12  and  DTU10  lead  to  lower  variance  for  six  missions  and 
GOT4.8  and  OSU12  for  four  missions,  while  other  models  lead  to  slightly  higher  variance  (differences  of 
0.05-0.42  cm).  These  results  suggest  that  all  seven  models  have  comparably  accurate  determination  of  the 
main  waves  when  averaged  over  the  global  ocean.  Certainly  there  remain  significant  differences  in  some 
local/shallow-water  regions  as  previously  shown  in  the  tide  gauges  tests  of  section  4,  but  these  differences 
are  for  the  most  part  evidently  smoothed  out  by  the  global  statistics  of  Table  1 1. 

7.  Elevation  Errors  of  Purely  Hydrodynamic  Models 

Although  this  paper  is  focused  mostly  on  highly  accurate  data-assimilative  global  tide  models,  it  seems  both 
timely  and  enlightening  to  use  the  test  machinery  we  have  in  place  to  examine  also  a  suite  of  numerical 
models  that  do  not  assimilate  any  data  at  all.  Such  models  are  far  less  accurate  than  those  that  assim¬ 
ilate  satellite  altimetry,  but  they  have  nonetheless  been  improving  significantly  over  the  past  decade. 
Improvements  in  purely  hydrodynamic  models  have  stemmed  from  higher  spatial  resolutions  and  better 


Table  1 1 .  Standard  Deviations  of  Residual  SSH  Anomaly  (in  cm)  Along  Satellite  Tracks  in  Deep  Oceans  (Depth 
>  1000  m)  After  Applying  Ocean  Tide  Model  Corrections3 


Model 

T/P 

Interlaced 

GFO 

Envisat 

Jason-1 

Interlaced 

Jason-1 

Jason-2 

Jason-1 

GM 

CryoSat-2 

GOT4.8 

10.07 

9.76 

10.30 

9.99 

10.12 

9.84 

10.06 

10.15 

10.21 

OSU12 

10.13 

9.60 

10.23 

9.97 

10.32 

9.76 

10.16 

10.34 

10.44 

DTU10 

10.07 

9.76 

10.31 

10.01 

10.13 

9.84 

10.07 

10.15 

10.22 

EOTlIa 

10.08 

9.75 

10.31 

9.98 

10.11 

9.84 

10.06 

10.16 

10.21 

HAM12 

10.07 

9.80 

10.32 

10.00 

10.12 

9.85 

10.06 

10.13 

10.19 

FES12 

10.10 

9.78 

10.32 

9.99 

10.13 

9.85 

10.07 

10.15 

10.21 

TPX08 

10.09 

9.80 

10.34 

10.02 

10.13 

9.87 

10.07 

10.14 

10.21 

Before 

32.65 

32.62 

33.00 

32.77 

32.62 

32.55 

32.53 

32.39 

32.92 

3The  bottom  line  is  the  respective  SSH  standard  deviation  before  ocean  tide  correction  was  applied. 


STAMMER  ETAL. 


©2014.  American  Geophysical  Union.  All  Rights  Reserved. 


266 


IgAGU  Reviews  of  Geophysics 


10.1 002/201 4RG000450 


Table  12.  Differences  (cm)  of  M2  Elevations  in  NSWC 
and  in  Purely  Hydrodynamic  Models  Versus  Test  Data3 


Deep  Ocean 

Shelf  Seas 

RMStg 

RMSalt 

RMStg 

RMSalt 

NSWC 

4.27 

4.41 

— 

17.4 

HIM 

8.75 

5.25 

33.7 

22.3 

OTIS-GN 

7.54 

6.76 

25.3 

18.6 

STORMTIDE 

8.33 

7.76 

48.2 

27.9 

OTIS-ERB 

5.63 

4.65 

23.6 

24.0 

STM-1  B 

12.69 

7.74 

30.5 

25.8 

HYCOM 

7.82 

7.00 

49.0 

26.2 

aRMSTG,  root-mean-square  differences  with  BPR 
and  tide  gauge  constants.  RMSALT,  root-mean-square 
differences  with  model  TPX08  over  the  entire 
ocean  domain. 


bathymetric  data  and  also  from  better  understand¬ 
ing  of  tidal  dissipation — in  particular  from  the  real¬ 
ization  that  open-ocean  conversion  of  barotropic 
tidal  energy  into  baroclinic  motion  is  a  significant 
energy  sink  [Egbert  and  Ray,  2000],  As  noted  in 
section  2.3,  all  but  one  of  the  six  forward  models 
examined  here  account  for  this  energy  sink  through 
a  parameterized  topographic  internal  wave  drag; 
Green  and  Nycander  [2013]  discuss  such  parameter- 
izations  in  some  detail.  Thus,  recent  improvements 
in  hydrodynamic  tide  models  represent  in  large 
part  a  better  understanding  of  the  ocean,  and  these 
advances  deserve  to  be  documented.  We  restrict  our 
attention  here  to  only  M2. 


Table  1 2  summarizes  comparisons  of  each  forward 
tide  model  with  the  same  sets  of  deep-ocean  bottom  pressure  recorders  and  shelf  water  tide  stations  used 
in  sections  4.1  and  4.2,  respectively.  As  anticipated,  all  the  tabulated  differences  for  the  forward  models  are 
much  larger,  by  roughly  a  factor  of  1 0,  than  corresponding  differences  for  the  modern  models  that  assimi¬ 
late  altimetry.  In  light  of  this,  for  testing  the  forward  models  we  may  treat  the  assimilation  models  as  "truth" 
and  thereby  obtain  a  more  complete  picture  of  the  modeling  errors,  one  not  confined  to  a  sparse  tide  gauge 
network.  Thus,  the  corresponding  RMS  differences  with  TPX08  have  been  added  to  Table  1 2.  The  tide  gauge 
and  TPX08  statistics  are  not  completely  consistent  in  their  ranking  of  models,  probably  owing  to  the  geo¬ 
graphic  limitations  of  the  gauge  network,  but  they  do  generally  agree  on  which  models  are  most  accurate 
and  which  are  less  accurate.  The  table  entries  suggest  that  the  forward  models  HIM  and  OTIS-ERB  are  the 
most  accurate  in  deep  water,  with  OTIS-GN  not  far  behind;  OTIS-GN  is  the  most  accurate  in  shallow  water, 
with  HIM,  STM-1  B,  and  OTIS-ERB  somewhat  less  accurate  and  the  others  trailing. 


It  is  interesting  to  observe  that  some  of  the  forward  models  are  now  starting  to  reach  accuracies  compa¬ 
rable  to  those  of  the  pre-Topex  era  assimilative  NSWC  model  of  Schwiderski  [1979],  which  from  the  late 
1 970s  to  early  1 990s  was  generally  considered  the  most  accurate  tidal  model  available.  The  RMS  differences 
between  NSWC  and  TPX08  are  4.41  cm  in  deep  water  and  1 7.4  cm  in  shallow  water,  which  is  only  slightly 
better  than  some  of  the  entries  in  Table  12.  For  models  that  lack  any  data  constraints  other  than  bathymetry, 
to  approach  the  NSWC  accuracies  does  mark  substantial  progress  in  modeling. 


Figure  1 2  shows  the  vector  differences  between  each  of  the  six  forward  models  and  TPX08.  To  put  these 
results  into  context,  the  vector  differences  between  the  altimeter-constrained  TPX08  and  GOT4.8  models, 
as  well  as  between  TPX08  and  NSWC,  are  included  in  the  figure.  As  was  to  be  expected,  the  TPX08-G0T4.8 
differences  are  almost  negligible  in  the  deep  ocean.  Consistent  with  Table  1 2,  the  NSWC  model  is  only 
slightly  more  accurate  than  some  of  the  forward  models.  According  to  Figure  12  (we  note  its  very 
nonlinear  color  bar),  all  six  forward  models  have  relatively  large  errors  in  the  North  Atlantic.  This  is  consistent 
with  the  large  M2  amplitudes  in  that  basin  and  the  Atlantic's  near-resonant  response  to  the  astronomi¬ 
cal  potential.  OTIS-GN  appears  the  most  accurate  there,  with  errors  exceeding  1 0  cm  only  in  the  very  high 
amplitude  areas  (e.g.,  Hudson  Strait). 

Three  models  do  not  include  the  ice-covered  Ross  and  Weddell  Seas,  and  all  three  are  seen  to  have  large 
errors  extending  throughout  much  of  the  Southern  Ocean;  HYCOM  in  particular  seems  sensitive  to  those 
missing  seas.  At  one  time  tidal  flexing  of  ice  shelves  was  conjectured  to  be  an  important  tidal  energy  sink, 
but  it  is  now  known  that  the  entire  Antarctic  coastline  accounts  for  less  than  about  50  GW  of  dissipation 
[Ray  and  Egbert,  1 997;  Egbert  and  Ray,  2001  ]  or  less  than  a  few  percent  of  the  global  total.  Nonetheless,  it 
seems  evident  that  incorrect  boundary  conditions  along  that  coastline  are  still  capable  of  inducing  model¬ 
ing  errors  into  far  lower  latitudes,  so  inclusion  of  these  Antarctic  seas  seems  critical  for  models  that  do  not 
assimilate  data. 


8.  Barotropic  Tidal  Currents 

Comparisons  of  modeled  tidal  currents  with  observations  are  much  less  common  than  comparisons  of 
modeled  tidal  elevations  with  observations,  especially  on  basin  and  global  scales.  Among  the  basin-  and 
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Figure  12.  Absolute  magnitudes  (cm)  of  vector  differences  in  M2  elevations  between  TPX08  and  purely  hydrodynamic 
("forward")  models  identified  by  the  labels  in  each  panel.  In  the  top  row  respective  differences  relative  to  GOT4.8  and 
NSWC  are  shown  for  reference. 


global-scale  studies  that  have  been  published,  three  representative  examples  are  Luyten  and  Stommel 
[1991],  who  compared  M2  currents  in  moored  current  meter  observations  with  those  in  the  Schwiderski 
[1979]  model;  Dushaw  et  al.  [1997],  who  compared  currents  from  several  tidal  constituents  in  TPXO2.0 
[Egbert  et  al.,  1 994]  to  currents  derived  from  acoustic  tomography  and  from  moored  current  meters;  and  Ray 
[2001],  who  deduced  tidal  currents  from  altimetrically  derived  tidal  elevations  and  compared  resulting  cur¬ 
rents  with  those  derived  from  models,  acoustic  tomography,  and  moored  current  meters.  All  three  of  these 
studies  focused  on  the  barotropic  component  of  tidal  currents.  Timko  et  al.  [2012,  201 3]  compared  fully 
three-dimensional  (barotropic  plus  baroclinic)  tidal  currents  in  global  eddy-resolving  HYCOM  with  those 
in  current  meters.  It  should  be  noted  that  most  barotropic  tidal  models  do  not  solve  for  currents  directly, 
so  much  as  transports,  from  which  currents  are  then  derived.  The  accuracy  of  model  currents  depends  on 
the  accuracy  of  local  seafloor  topography  employed  for  the  model.  As  noted  above,  while  some  of  the  tidal 
models  are  constrained  by  observations  of  tidal  surface  elevation,  none  of  the  tidal  models  are  constrained 
by  measurements  of  tidal  currents.  While  tidal  currents  are  inherent  to  hydrodynamic  models,  some  of  the 
empirical  models  consider  elevations  only,  with  currents  computed  after  the  fact. 
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In  this  section,  M2  barotropic  tidal  currents  from  both  data-assimilative  models  and  forward  models  are 
compared  against  those  derived  from  historical  moored  current  meter  records  (section  8.1)  and  from  acous¬ 
tic  tomography  (section  8.2).  The  data-assimilative  models  used  here  are  the  TPX08,  GOT4.8,  HAM1 2,  FES1 2, 
and  NSWC  models.  The  forward  models  used  here  are  the  barotropic  HIM,  OTIS-ERB,  OTIS-GN,  and  STM-1  B 
models,  and  the  baroclinic  HYCOM  and  STORMTIDE  models. 

8.1.  Tidal  Currents  From  Moored  Current  Meters 

The  moored  current  meter  records  utilized  here  represent  a  subset  of  those  employed  by  Scott  et  at.  [2010] 
and  Timko  et  at.  [2012,  2013].  Because  the  present  study  focuses  on  barotropic  tidal  currents,  we  chose 
only  moorings  for  which  a  meaningful  vertical  average  might  be  computed.  We  required  sampling  at  a 
minimum  of  once  per  hour,  for  65  days,  and  at  a  minimum  of  three  vertical  levels,  with  at  least  one  above 
800  m  and  one  below  1 950  m.  After  excluding  moorings  with  especially  obvious  measurement  problems 
at  small  velocities  (more  below),  or  other  apparent  problems,  a  total  of  56  moorings  (Figure  1 3a)  met 
these  requirements. 

The  geographic  sampling  of  tidal  currents  by  moorings  is  very  inhomogeneous  given  that  mooring  loca¬ 
tions  tend  to  focus  on  western  boundary  current  regions,  leaving  large  regions,  such  as  the  entire  South 
Pacific,  undersampled.  Vertical  coverage  is  also  limited.  Of  the  56  moorings  used  here,  24  have  only  three 
current  meters,  1 9  have  only  four,  8  have  only  five,  and  only  5  have  six  or  more.  The  sparse  vertical  sampling 
of  the  moorings  makes  accurate  estimation  of  the  barotropic  mode  difficult,  with  considerable  contamina¬ 
tion  from  baroclinic  modes,  as  noted  in  earlier  studies  [e.g.,  Luyten  and  Stommel,  1991;  Dushaw  et  at.,  1997]. 
Indeed,  some  earlier  studies  [e.g.,  Luyten  and  Stommel,  1991],  instead  of  estimating  the  barotropic  mode 
directly,  took  the  currents  in  the  bottommost  current  meter  of  a  mooring  to  represent  the  barotropic  mode. 
The  HYCOM  simulations,  with  relatively  high  vertical  resolution,  will  be  used  to  estimate  the  sampling  error 
in  M2  barotropic  tidal  currents  computed  from  vertically  sparse  moored  current  meter  records. 

We  separate  the  barotropic  and  baroclinic  components  of  tidal  currents  at  the  mooring  locations  using  stan¬ 
dard  linear  dynamical  modes  Fm(z),  where  m  is  a  vertical  mode  index  and  z  denotes  depth.  The  modes  are 
obtained  from  the  internal  wave  Sturm-Liouville  problem  [e.g.,  Gill,  1982;  Wunsch,  1975],  written  in  terms 
of  the  Brunt-Vaisala  buoyancy  frequency  N(z )  [Gill,  1 982],  and  an  internal  wave  frequency  taken  to  be  that 
of  M2.  We  take  the  baroclinic  modes  to  begin  with  an  index  of  m  =  2  (for  the  first  baroclinic  mode);  the 
barotropic  mode  F^z)  =  1  at  all  depths.  For  each  mooring  location,  A/(z)  is  computed  from  the  nearest 
grid  point  of  the  climatological  2009  World  Ocean  Atlas  [Antonov  et  at.,  201 0;  Locarnini  et  at.,  201 0].  For 
HYCOM,  A/(z)  is  computed  from  the  model  output  itself.  The  vertical  modes  are  obtained  with  a  shooting 
method  code  (Glenn  Flierl,  personal  communication,  1 995),  which  yields  results  that  are  orthonormal  to 
high  precision. 

We  can  write  the  zonal  component  u  of  either  the  moored  current  meter  velocities  or  the  HYCOM  velocities 
sampled  at  the  moored  current  meter  vertical  levels  (referred  to  as  "SAMPLED  HYCOM"  in  what  follows),  as 

K  L 

U(Zl  f)  =  Z  Z  ref)  C0S  [®c(f  “  fref)  +  Xc« ref)  +  vc(fref)  “  C]  +  6’  (4) 

m= 1  c=  1 

where  f  is  time,  tref  is  a  reference  time  (i.e.,  the  beginning  of  the  first  year  of  a  mooring  record),  m  is  an  index 
of  vertical  mode  number,  K  is  the  total  number  of  vertical  modes  employed  including  the  barotropic  mode, 
c  is  an  index  of  tidal  constituent  (M2,  S2,  K,,  etc.),  L  is  the  total  number  of  constituents  used,  fc(tref )  and  vc(fref) 
are  respectively  the  constituent-dependent  nodal  modulations  of  amplitude  and  phase  referenced  to  tref,  a>c 
is  the  frequency  of  constituent  c,  jc(fref )  is  the  constituent-dependent  astronomical  argument  referenced  to 
tref  [e.g.,  Pugh,  1 987],  and  e  is  a  residual. 

The  mode-  and  constituent-dependent  amplitudes  Acm  and  phase  lags  <pcm  are  solved  for  with  a  linear  least 
squares  fit.  In  effect,  we  have  extended  standard  tidal  harmonic  analysis  to  solve  for  vertical  mode  structure 
over  all  the  current  meters  in  a  mooring.  A  similar  problem  is  solved  for  the  meridional  velocity  compo¬ 
nent  v.  We  solved  for  the  five  largest  tidal  constituents  ( L  =  5;  M2,  S2,  N2,  K, ,  and  O, )  but  found  that  the  M2 
barotropic  mode  estimated  from  the  above  technique  was  not  very  sensitive  to  the  inclusion  of  other  con¬ 
stituents.  We  found  a  slightly  closer  agreement  between  the  current  meter  and  TPX08  estimates  of  the  M2 
barotropic  tidal  currents  with  two  vertical  modes  in  the  mooring  estimation  than  with  three.  In  what  follows 
we  display  results  from  a  fit  using  only  two  vertical  modes  (barotropic  and  first  baroclinic). 
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inferred  from  mooring  data,  (c)  M2  tidal  ellipses  inferred  from  HAM12  at  the  mooring  positions.  In  Figures  13b  and  13c 
only  39  locations  are  shown  to  avoid  clutter  due  to  moorings  located  close  to  each  other. 


Figure  1 3b  shows  the  M2  barotropic  tidal  current  ellipses  computed  from  the  moorings  via  the  procedure 
above.  Figure  1 3c  shows  the  M2  barotropic  tidal  current  ellipses  computed  from  the  altimeter-constrained 
model  FIAM12.  Visual  comparison  of  Figures  13b  and  13c  indicates  that  there  are  many  locations  where  the 
mooring  and  FIAM1 2  estimates  lie  close  to  one  another  but  also  many  locations  where  they  do  not.  There 
are  some  locations  where  the  currents  are  of  large  magnitude  and  are  quite  different  in  the  two  estimates. 
This  suggests  that  relatively  small  numbers  of  outliers  may  control  RMS  statistics,  in  like  manner  to  earlier 
discussions  of  model/tide  gauge  elevation  differences. 
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Figure  14.  Scatterplots  of  (a,  b)  amplitudes  and  (c,  d)  phase  lags  of  M2  barotropic  tidal  zonal  velocity  component  u 
at  moored  current  meter  locations  shown  in  Figure  13,  estimated  from  five  different  data-constrained  barotropic  tide 
models,  four  different  forward  barotropic  tide  models,  and  two  different  forward  multilayer  tide  models  forced  by  atmo¬ 
spheric  fields  in  addition  to  the  astronomical  tidal  potential,  all  plotted  versus  the  mooring  results  (Figures  14a  and  14c) 
and  data-constrained  model  TPX08  (Figures  14b  and  14d).  Values  from  data-constrained  and  forward  barotropic  models 
denoted  by  crosses  and  circles,  respectively.  Values  from  multilayer  forward  models  denoted  by  squares. 

Figure  14  shows  the  amplitudes  and  phase  lags  of  the  M2  barotropic  zonal  velocity  component  (u)  from  the 
various  assimilative  and  forward  tide  models,  plotted  against  the  mooring  estimates  (Figures  14a  and  14c) 
and  against  the  TPX08  altimeter-constrained  model  (Figures  14b  and  14d).  For  both  amplitude  and  phase 
lag  there  is  more  scatter  when  models  are  plotted  against  mooring  values  than  against  TPX08  values.  This 
is  consistent  with  previous  work  [e.g.,  Dushaw  et  a!.,  1997;  Ray,  2001],  which  concluded  that  tidal  currents  in 
different  models  often  agree  more  closely  with  each  other  than  with  estimates  made  from  moorings.  The 
scatter  in  amplitude  is  somewhat  larger  than  in-phase  lag.  The  equivalent  scatterplot  for  the  meridional 
velocity  component  v  (Figure  1 5)  behaves  similarly  but  shows  somewhat  less  scatter  than  the  u  plots. 

Because  barotropic  velocity  consists  of  two  (zonal  and  meridional)  scalar  components,  we  define  the  RMS 
velocity  difference  D  between  two  different  estimates  as  a  quadrature  sum  of  the  RMS  differences  for  each 
separate  component;  thus,  the  RMS  velocity  difference  D  is  a  two-component  velocity  equivalent  of  the  RMS 
elevation  difference  which  was  used  extensively  earlier  in  the  paper.  As  with  elevations  (see,  e.g.,  Table  3) 
we  compute  a  signal  strength  S  alongside  the  RMS  difference  D;  again,  as  with  D,  the  velocity  signal  is  a 
two-component  equivalent  to  the  elevation  signal  used  earlier. 

The  RMS  differences  are  computed  against  two  standards,  the  moored  current  meter  archive  (CMA)  and 
TPX08.  All  resulting  signals  S  and  RMS  differences  DCMA  and  DTPX0  are  listed  in  Table  1 3.  While  the  historic 
Schwiderski  [1979]  NSWC  model  has  by  far  the  largest  signal  S,  values  from  all  modern  models  lie  relatively 
close  together,  ranging  from  1.53  cm/s  (OTIS-ERB)  to  1.89  cm/s  (STORMTIDE),  a  range  of  about  24%.  Assim¬ 
ilative  models  lie  in  a  tighter  range,  from  1 .57  cm/s  (GOT4.8)  to  1.67  cm/s  (FIAM12) — a  range  of  about  7%. 
Table  1 3  confirms  that,  consistent  with  discussions  in  previous  studies,  the  signal  in  the  moored  current 
meters  is  slightly  weaker  than  the  signals  in  any  of  the  models.  Luther  et  at.  [1991]  and  other  sources  discuss 
instrumental  reasons  (such  as  stalled  rotors)  that  contribute  to  this  tendency.  Our  own  visual  inspection  of 
the  bottommost  current  meter  records  in  the  CMA  indicates  that  in  approximately  one  third  of  the  mooring 
locations  examined,  there  are  clear  measurement  problems  at  small  velocities,  consistent  with  stalled  rotors. 
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Mooring  v  amplitude  (cm/s)  TPX08  v  amplitude  (cm/s) 


ra 


■§.  -100“ — - - — 1 1 

>  0  200  400 

Mooring  v  phase  lag  (degrees) 


o  w 

C  <D 
d>  0) 


400 
300 
200 
100 
0 

-100 

TPX08  v  phase  lag  (degrees) 


(d) 


200  400 


X 

TPX08 

X 

GOT4.8 

X 

HAMTIDE 

X 

FES 

X 

NSWC 

O 

STM-1  B 

O 

OTIS-GN 

O 

HIM 

OTIS-ERB 

□ 

STORMTIDE 

□ 

HYCOM 

Figure  15.  Same  as  Figure  14  but  for  meridional  velocity  component  v. 


The  RMS  difference  of  SAMPLED  HYCOM  versus  fully  vertically  sampled  HYCOM  is  0.99  cm/s.  This  num¬ 
ber  roughly  represents  an  error  in  barotropic  mode  estimation  from  sparse  vertical  sampling.  We  therefore 
cannot  reasonably  expect  any  of  the  models  to  agree  with  the  CMA  estimates  at  a  level  much  lower  than 
0.99  cm/s.  Indeed,  the  DCMA  values  are  of  this  order,  generally  a  little  less  for  modern  assimilative  models, 
a  little  more  for  modern  forward  models,  and  much  more  for  the  NSWC  and  SAMPLED  HYCOM  results.  The 


Table  13.  Signals  S  and  RMS  Differences  D  of  M2  Tidal  Currents  (See  Text  for 
Definition),  Averaged  Over  the  Mooring  Locations  Shown  in  Figure  1 3a 


Model  or  Data  Set 

Signal  (cm/s) 

Dcma  (cm/s) 

f)TPX08  (cm/s) 

Mooring  data 

1.45 

- 

0.83  (0.39) 

Modern  data-constrained  models 

TPX08 

1.63 

0.83  (0.39) 

GOT4.8 

1.57 

0.87  (0.38) 

0.32  (0.047) 

HAM12 

1.67 

0.90  (0.38) 

0.24  (0.049) 

FES12 

1.65 

0.83  (0.40) 

0.16  (0.086) 

Historical  Models 

NSWC 

2.24 

1.58  (0.38) 

1.50  (0.23) 

Purely  hydrodynamic  models 

HIM 

1.59 

1.07  (0.56) 

0.65  (0.17) 

OTIS-GN 

1.56 

0.98  (0.57) 

0.47  (0.22) 

STORMTIDE 

1.89 

1.18  (0.55) 

0.62  (0.30) 

OTIS-ERB 

1.53 

1.03  (0.51) 

0.52  (0.20) 

STM-1  B 

1.72 

1.17  (0.64) 

0.73  (0.54) 

HYCOM 

1.86 

1.06  (0.67) 

0.56  (0.38) 

SAMPLED  HYCOM 

1.87 

1.52  (1.07) 

1 .05  (0.84) 

aRMS  differences  DCMA  and  Dypxos  are  computed  against  the  moored  current 
meter  archive  (CMA)  and  against  TPX08,  respectively.  The  RMS  difference  of  SAM¬ 
PLED  HYCOM  versus  HYCOM  is  0.99  cm/s.  Numbers  in  parentheses  are  median 
D  values. 
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Figure  16.  M2  barotropic  tidal  current  ellipses  for  (top)  the  central  North  Pacific  and  (bottom)  the  western  subtropical 
Atlantic  derived  from  the  TPX08  model.  The  RTE87  (1987),  AMODE  (1991),  and  HOME  (2001)  tomography  arrays  are 
indicated  in  red.  Barotropic  currents  were  estimated  along  these  paths  using  reciprocal  acoustic  transmissions. 
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Figure  17.  (a-f)  Polar  plot  of  M2  tidal  current  amplitude  and  phase  derived  from  various  tide  models  (black: 
non-data-assimilating,  red:  assimilating,  cyan:  Schwiderski)  as  averages  along  the  tomography  paths  indicated  in  each 
panel.  Tomographic  estimates  are  provided  as  green  error  bars  (see  Tables  14  and  15).  Results  for  the  Schwiderski  model 
shown  in  Figures  17c  and  17d  were  reported  by  Dushaw  et  at.  [1995]  and  recomputed  in  this  paper,  respectively. 


slightly  lower  DCMA  values  for  modern  assimilative  versus  forward  tide  models  suggest  that  the  barotropic 
tidal  current  estimates  made  from  moorings  retain  some  minimal  power  to  discriminate  between  the 
qualities  of  the  two  different  types  of  models. 

From  the  above  it  can  be  expected  that  the  largest  DTPX08  values  in  Table  1 3  can  be  found  for  NSWC,  SAM¬ 
PLED  HYCOM,  and  the  current  meter  archive.  The  DTPX08  values  for  all  other  models  range  from  0.16  cm/s 
(FES1 2)  to  0.73  cm/s  (STM-1  B),  smaller  than  the  RMS  differences  found  relative  to  the  current  meter  archive. 
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Table  1 4.  RTE87  Array:  Tidal  Current-Velocity  Comparison,  Amplitudes,  and  Phases3 


Path  3 

->  1 

Path  2 

1 

Path3 

->  2 

Amplitude  (cm  s  ') 

Phase  Lag  (deg) 

Amplitude  (cms-1) 

Phase  Lag  (deg) 

Amplitude  (cm  s  ’) 

Phase  Lag  (deg) 

AT 

TP 

AT 

TP 

AT 

TP 

AT 

TP 

AT 

TP 

AT 

TP 

1.31  ±0.03 

1.24 

223  ±1 

221 

1.54  ±0.04 

1.57 

176  ±2 

173 

0.95  ±0.07 

0.97 

131  ±4 

138 

S2 

0.53 

0.56 

272  ±4 

272 

0.56 

0.63 

207  ±6 

202 

0.33 

0.53 

1 70  ±14 

163 

n2 

0.14 

0.16 

191  ±13 

191 

0.35 

0.29 

169  ±8 

163 

0.16 

0.16 

111  ±23 

146 

k2 

0.12 

0.15 

268  ±  1 3 

269 

0.15 

0.18 

1 84  ±16 

196 

0.16 

0.15 

143  ±20 

158 

O, 

0.43 

0.45 

101  ±4 

105 

0.25 

0.26 

207  ±9 

204 

0.13 

0.30 

255  ±23 

253 

Ki 

0.75 

0.70 

128  ±  2 

125 

0.41 

0.39 

216  ±  7 

210 

0.31 

0.42 

306  ±  1 3 

273 

Pi 

0.18 

0.22 

1 32  ±11 

123 

0.14 

0.12 

229  ±21 

208 

0.07 

0.14 

171  ±61 

267 

Qi 

0.10 

0.08 

64±  16 

93 

0.06 

0.05 

1 88  ±  34 

202 

0.11 

0.06 

21 7  ±  28 

249 

3  AT  and  TP  stand  for  acoustic  tomography  and  TPX08,  respectively.  Figure  16  shows  the  locations  of  the  RTE87  moorings. 


The  RMS  differences  of  modern  assimilative  models  versus  TPX08  are  significantly  smaller  than  the  RMS 
differences  of  forward  models  versus  TPX08.  The  smaller  DTPX08  values  seen  in  modern  models,  both  for¬ 
ward  and  assimilative,  and  the  larger  difference  of  the  older  NSWC  model  versus  TPX08  indicate  a  likely 
improvement,  or  at  least  convergence,  overtime  of  currents  in  modern  tide  models  with  the  advent  of  satel¬ 
lite  altimetry  (for  the  assimilative  models)  and  topographic  wave  drag  schemes  (for  most  of  the  forward  tide 
models,  STORMTIDE  being  the  exception).  It  should  be  noted,  however,  that  even  the  relatively  small  RMS 
differences  of  the  modern  assimilative  tide  models  measured  against  TPX08  represent  relative  errors  (mea¬ 
sured  with  respect  to  the  signal  S)  of  about  10-20%,  much  larger  than,  for  instance,  the  2%  relative  errors  in 
M2  elevations  in  modern  assimilative  models  measured  with  respect  to  tide  gauges  (e.g.,  Table  3).  The  larger 
relative  errors  in  tidal  current  estimates  compared  to  tidal  elevation  estimates  are  another  indicator  of  the 
greater  difficulty  in  estimating  and  validating  tidal  currents  as  compared  to  tidal  elevations. 

Finally,  we  note  that,  as  in  the  elevation  comparisons  discussed  earlier  in  the  paper,  the  median  differences 
(given  inside  parentheses  in  Table  1 3)  are  lower  than  the  differences  defined  as  RMS  values.  In  the  case 
of  the  NSWC  model,  the  median  differences  are  much  lower  than  the  RMS  difference;  the  large  change 
seen  in  the  two  metrics  is  explained  by  very  large  differences  at  just  two  of  the  56  locations.  With  this 
major  exception,  the  general  patterns  seen  in  the  RMS  differences  are  also  seen  in  the  median  values;  the 
data-assimilative  models  perform  better  than  the  forward  models,  and  the  models  lie  closer  to  each  other 
than  to  the  mooring  estimates. 

8.2.  Tidal  Currents  From  Acoustic  Tomography 

Barotropic  tidal  currents  are  also  measured  by  reciprocal  acoustic  tomography,  which  employs  the  travel 
times  of  acoustic  rays  that  cycle  throughout  the  water  column  and  traverse  O  (500  km)  ranges  [Dushaw  et 
al.,  1 994].  The  measurement  is  inherently  integral  over  both  depth  and  range  and  is  not  significantly  influ¬ 
enced  by  baroclinic  tidal  currents.  Previous  comparisons  of  tidal  current  harmonic  constants  to  tidal  models 
have  shown  that  tomography  provides  precise  measurements  of  at  least  the  first  eight  tidal  constituents; 
M2,  S2,  N2,  K2,  K,,  O,,  and  even  P,  and  Q,  [Dushaw  et  al.,  1995, 1997].  Our  comparisons  here  will  employ  tidal 
currents  determined  from  three  experiments  in  the  North  Pacific  and  North  Atlantic  (Figure  16):  (1)  In  the 
central  North  Pacific,  tidal  currents  were  measured  along  the  three  legs  of  the  tomography  triangle  of  the 
Reciprocal  Tomography  Experiment  during  summer  1 987  (RTE87)  [ Dushaw  et  al.,  1 995].  (2)  The  tomography 
arrays  deployed  in  the  vicinity  of  the  Fiawaiian  Ridge  in  2001  as  part  of  the  Flawaiian  Ocean  Mixing  Experi¬ 
ment  (FIOME)  [ Rudnick  et  al.,  2003;  Dushaw  et  al.,  2011]  were  intended  to  measure  barotropic  tidal  currents 
with  concurrent  pressure  measurements.  The  arrays  were  also  designed  to  measure  radiation  of  internal 
tides  from  the  Flawaiian  Ridge  [Dushaw  et  al.,  201 1];  hence,  the  longest  path  of  the  diamond-shaped  arrays 
were  parallel  to  the  Ridge.  (3)  In  the  western  North  Atlantic,  tidal  currents  were  measured  along  each  of  the 
15  legs  of  the  Acoustic  Mid-Ocean  Dynamics  Experiment  (AMODE)  tomography  array  in  1991-1992  [ Dushaw 
et  al.,  1 997].  The  AMODE  array  had  path  lengths  of  350-660  km,  1 50-250  day  record  lengths,  and  mea¬ 
surements  obtained  every  3  h  on  every  fourth  day.  The  tomography  and  model  estimates  for  tidal  currents 
along  several  tomography  paths  are  compared  in  Figure  17  and  Tables  14  and  15.  The  figure  shows  only  a 
selection  of  comparisons,  and  the  tables  include  only  a  comparison  relative  to  just  the  TPX08  model  due  to 
space  limitations. 
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Table  1 5.  HOME  Arrays:  Tidal  Current-Velocity  Comparisons,  Amplitudes,  and  Phases3 
South  Path  1  -»  2  South  Path  1  ->  3 


Amplitude  (cm  s  ’)  Phase  Lag  (deg)  Amplitude  (cm  s  ')  Phase  Lag  (deg) 


AT 

TP 

AT 

TP 

AT 

TP 

AT 

TP 

m2 

1 .09  ±  0.03 

1.20 

152  ±2 

149 

1.00  ±0.02 

1.21 

207  ±1 

208 

S2 

0.42 

0.51 

182  ±4 

167 

0.47 

0.56 

229  ±3 

230 

n2 

0.20 

0.22 

163  ±  9 

146 

0.21 

0.23 

204  ±6 

207 

k2 

0.12 

0.14 

1 82  ±15 

161 

0.16 

0.15 

212  ±  7 

224 

O, 

0.14 

0.19 

236  ±13 

245 

0.23 

0.27 

182  ±5 

178 

K, 

0.17 

0.17 

239  ±  10 

258 

0.46 

0.49 

173  ±3 

178 

Pi 

0.07 

0.06 

31 7  ±26 

257 

0.14 

0.15 

171  ±8 

180 

Qi 

0.03 

0.04 

227  ±59 

240 

0.04 

0.04 

275  ±28 

177 

South  Path  1  -»  4 

South  Path  2  ->  3 

Amplitude  (cm  s  ') 

Phase  Lag  (deg) 

Amplitude  (cm  s  ') 

Phase  Lag  (deg) 

AT 

TP 

AT 

TP 

AT 

TP 

AT 

TP 

m2 

1.15  ±0.03 

1.33 

263  ±1 

264 

1.1 8  ±0.04 

1.29 

285  ±2 

278 

S2 

0.61 

0.61 

280  ±3 

280 

0.59 

0.60 

292  ±4 

294 

n2 

0.27 

0.27 

264  ±6 

260 

0.25 

0.25 

276  ±9 

273 

k2 

0.16 

0.17 

276±  10 

276 

0.19 

0.17 

31 8  ±  1 3 

290 

O, 

0.33 

0.29 

133  ±5 

140 

0.23 

0.26 

1 23  ±10 

127 

K, 

0.49 

0.54 

166±  3 

161 

0.43 

0.46 

153  ±5 

154 

Pi 

0.12 

0.16 

1 43  ±14 

160 

0.20 

0.14 

138  ±12 

153 

Q, 

0.04 

0.05 

99  ±46 

124 

0.09 

0.05 

14±27 

109 

South  Path  2  -»  4 

South  Path  3  -►  4 

Amplitude  (cm  s  ') 

Phase  Lag  (deg) 

Amplitude  (cm  s  ') 

Phase  Lag  (deg) 

AT 

TP 

AT 

TP 

AT 

TP 

AT 

TP 

m2 

1.22  ±0.02 

1.32 

298  ±1 

294 

1.33  ±0.03 

1.36 

311  ±1 

309 

S2 

0.55 

0.58 

311  ±2 

310 

0.44 

0.57 

318±4 

327 

n2 

0.23 

0.25 

286  ±5 

289 

0.25 

0.26 

311  ±8 

303 

k2 

0.14 

0.16 

308  ±9 

305 

0.13 

0.16 

31 7  ±  15 

322 

Oi 

0.21 

0.23 

1 1 3  ±  6 

110 

0.16 

0.20 

82  ±  12 

87 

Ki 

0.30 

0.35 

1 1 3  ±  4 

110 

0.13 

0.21 

89  ±  14 

124 

Pi 

0.11 

0.10 

1 57  ±11 

142 

0.02 

0.07 

9  ±  1 08 

120 

Q, 

0.04 

0.05 

1 10  ±  32 

93 

0.09 

0.04 

33  ±  21 

76 

North  Path  1  ->  3 

Amplitude  (cm  s  ’) 

Phase  Lag  (deg) 

AT 

TP 

AT 

TP 

m2 

1.55  ±0.02 

1.67 

10±1 

12 

S2 

0.73 

0.76 

38  ±2 

40 

n2 

0.30 

0.30 

12  ±4 

10 

k2 

0.25 

0.21 

28  ±5 

35 

Oi 

0.21 

0.25 

354  ±5 

11 

Ki 

0.32 

0.44 

350  ±4 

10 

vP, 

0.16 

0.13 

350  ±7 

10 

Qi 

0.03 

0.04 

13  ±40 

7 

3 AT  and  TP  stand  for  acoustic  tomography  and  TPX08,  respectively.  Figure  16  shows 
the  locations  of  moorings. 


In  the  North  Atlantic,  the  harmonic  constants  derived  from  the  AMODE  paths  along  the  major  (Path  2,  5) 
and  minor  axes  (Path  1,4)  of  the  tidal  ellipses  indicate  that  the  comparisons  between  measured  and  model 
tidal  harmonic  constants  are  in  general  agreement  (Figures  17a  and  1 7b).  Estimates  from  three  of  the  four 
data-assimilating  models  are  tightly  clustered  (the  outlier  is  FES1 2),  with  variations  in  amplitude  and  phase 
of  only  O  (0.05  cm/s)  and  a  few  degrees.  Not  surprisingly,  models  not  constrained  by  data  gave  estimates 
with  greater  scatter  but  were  still  in  agreement  with  the  other  models  and  tomography  by  O  (0.1  cm/s)  in 
amplitude  and  0(10°)  in  phase.  With  a  measurement  uncertainty  of  0.02  cm/s  and  a  difference  of  about 
0.2  cm/s,  the  tomography  amplitude  along  the  major  axis  disagrees  with  the  assimilation  models  by  about 
1 0  standard  deviations.  Except  for  the  FES1 2  model,  the  M2  tidal  phases  measured  on  all  the  tomography 
paths  agree  with  those  from  the  assimilation  models  within  uncertainty. 
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In  the  North  Pacific,  harmonic  constants  for  tidal  currents  from  data-assimilating  models  are  also  tightly 
clustered,  while  harmonic  constants  from  non-data-assimilating  models  are  noticeably  more  scattered 
(Figures  1 7c-1 7f).  In  the  central  North  Pacific,  RTE87  tomography  and  data-assimilating  models  agree  within 
measurement  uncertainty  (0.03  cm/s,  1  °).  Near  Hawaii  (Figures  17e  and  17f),  however,  systematic  differ¬ 
ences  between  measured  and  modeled  harmonic  constants  (0.25  cm/s  in  amplitude  and  a  few  degrees  in 
phase)  are  significant.  As  observed  in  the  Atlantic,  differences  occur  primarily  in  amplitude  and  not  phase. 
We  can  offer  no  explanation  for  the  differences  between  the  M2  harmonic  constants  estimated  from  tidal 
models  and  tomography  at  this  time. 

One  way  to  test  the  self-consistency  of  harmonic  constants  and  their  uncertainties  obtained  from  acous¬ 
tic  paths  of  a  triangular  tomographic  array  is  to  estimate  tidal  vorticity  [Dushaw  et  al.,  1 997].  By  Stokes' 
theorem,  integration  of  current  around  a  closed  loop  is  equal  to  the  areal-average  relative  vorticity.  For 
barotropic  tidal  currents,  this  relative  vorticity  is  small,  since  barotropic  currents  are  predominantly  irrota- 
tional.  The  small  vorticity  of  the  tides  is  induced  primarily  by  tidal  elevation  stretching  the  vortex  lines,  with 
lesser  contributions  through  meridional  flow  (the  "beta-v"  term)  and  flow  over  varying  topography.  The 
AMODE  tomographic  array  offers  1 0  triangles  from  which  vorticity  was  computed  [ Dushaw  et  al.,  1 997].  The 
relative  vorticities  computed  in  this  way  had  amplitudes  larger  than  those  obtained  from  the  models  by  a 
factor  of  about  3-4  and  had  fairly  stable  phases  that  led  most  models  by  10-50°.  Uncertainties  were  large 
at  about  20%.  This  calculation  indicates  that  the  harmonic  constants  estimated  on  the  various  tomography 
paths  are  at  least  self-consistent  with  the  formal  estimates  for  their  uncertainties.  The  estimated  relative  vor¬ 
ticity  was  larger  than  expected,  however,  suggesting  that  the  systematic  differences  between  tomography 
and  tidal  models  may  be  caused  by  larger-than-expected  rotational  components  to  the  tidal  currents. 

9.  Summary  and  Outlook 

As  demonstrated  in  this  study,  the  accuracy  of  tide  models  has  improved  substantially  over  the  last  20  years. 
As  an  example,  we  find  an  RMS  difference  of  M2  tide  models  relative  to  pelagic  bottom  pressure  stations  of 
0.51  cm  which  has  to  be  compared  with  1 .64  cm  reported  by  Shum97.  However,  about  half  of  this  apparent 
model  improvement  actually  stems  from  improvements  in  the  test  data,  highlighting  the  fact  that  even 
today  a  substantial  effort  needs  to  be  spent  on  the  observational  in  situ  tidal  data  bases.  (Our  equivalent  test 
of  the  best  Shum97  model  gives  0.92  cm,  but  this  is  still  roughly  twice  the  RMS  of  the  best  modern  model.) 
Meanwhile,  improvements  of  global  models  in  shallow  water  have  been  very  impressive,  with  M2  differences 
falling  from  23.0  cm  to  3.5  cm  (Table  4,  bottom).  In  some  cases  we  were  able  to  identify,  and  sometimes 
correct,  errors  in  historic  tide  gauge  data  by  observing  when  model-gauge  differences  far  exceed  the  spread 
in  our  ensemble  of  models. 

According  to  our  tide  gauge  and  bottom  pressure  tests,  the  RMS  accuracy  of  the  best  altimeter-constrained 
tide  models  is  found  for  the  M2  constituent  to  be  0.51, 3.5,  and  4.5  cm  for  pelagic,  shelf,  and  coastal  condi¬ 
tions,  respectively.  For  the  combined  eight  major  tidal  constituents,  the  root  sum  of  squares  is  0.89,  5.1,  and 
6.5  cm  for  pelagic,  shelf,  and  coastal  conditions,  respectively.  Note,  however,  that  RSS  and  RMS  values  for 
shallow  regions  may  be  misleadingly  pessimistic  because  they  can  be  skewed  high  by  just  one  or  two  diffi¬ 
cult  locations  where  the  models  are  especially  poor.  This  was  emphasized  by  Figure  3,  which  shows  that  the 
errors  at  a  majority  of  shelf  locations  are  less  than  2.5  cm  in  M2  but  that  there  are  remaining  locations  where 
errors  are  much  larger  by  many  decimeters  or  more.  Some  of  these  error-prone  regions  are  highlighted  in 
Figure  1  but  not  always — for  example,  Zanzibar,  Tanzania,  was  identified  as  a  poor  station  for  some  mod¬ 
els  even  though  Figure  1  suggests  that  the  east  coast  of  Africa  is  not  especially  problematic.  Thus,  while  our 
overall  RMS  values  are  likely  a  good  assessment  of  model  quality  in  the  deep  ocean,  no  single  value  is  repre¬ 
sentative  of  the  errors  in  shallow  water;  the  models  may  be  fairly  accurate  or  may  be  very  poor,  depending 
on  location. 

The  long-wavelength  components  of  models  tested  by  analyzing  satellite  laser  ranging  measurements 
and  GRACE  intersatellite  ranging  measurements  confirm  that  none  of  the  tested  tide  models  substantially 
stands  out  in  overall  accuracy.  Several  of  the  modern  models  are  of  comparable  quality  for  use  in  precise 
orbit  determination.  The  GRACE  data  show  that  all  data-constrained  models  are  still  imperfect  on  basin  and 
subbasin  scales  and  that  no  single  model  outperforms  all  other  models,  especially  in  high  latitudes  where 
insufficient  high-quality  tide  information  exists.  Of  the  tested  models,  GOT4.8,  HAMTIDE,  and  TPX08  tend  to 
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have  the  smallest  RSS  differences  relative  to  in  situ  observations  for  the  deep  ocean.  For  shelf  regions,  FES1 2 
and  TPX08  appear  to  perform  best.  Over  the  Arctic,  DTU10  shows  the  lowest  RSS  values. 

Assessing  model  accuracies  of  tidal  currents  remains  problematic  owing  to  the  relative  sparseness  of  in  situ 
velocity  measurements  and  the  difficulty  of  isolating  the  barotropic  mode  in  moored  current  meter  mea¬ 
surements.  Relative  to  moored  current  meter  data,  reciprocal  acoustic  tomography  leads  to  estimates  of 
harmonic  constants  that  have  much  lower  reported  uncertainties,  but  unfortunately  such  measurements 
have  been  made  in  only  a  few  locations.  Nevertheless,  both  current  meter  and  tomography  comparisons 
show  that  our  tested  assimilation  models  consistently  offer  more  accurate  estimates  for  currents  than  do 
unconstrained  models.  We  also  observed  small  systematic  differences  between  all  M2  models  and  the 
tomography  data  in  the  western  North  Atlantic  and  near  Flawaii,  with  currents  from  models  tending  to  be 
larger  by  10-20%.  No  significant  discrepancies  were  found  between  modeled  and  observed  phases.  An 
explanation  for  the  amplitude  differences  awaits  future  studies.  Increasing  the  number  of  in  situ  high  verti¬ 
cal  resolution  current  measurement  records,  which  presently  constitute  a  small  data  set,  would  improve  our 
knowledge  of  barotropic  (and  baroclinic)  tidal  velocity  fields. 

While  global  ocean  tide  models  have  reached  an  impressive  level  of  accuracy,  there  remain  outstanding 
challenges  that  need  to  be  overcome,  especially  in  shallow-water  environments  and  in  the  extreme  and 
high-latitude  environments.  The  principal  limitations  of  modern  tide  models  in  polar  regions  arise  from  a 
combination  of  poor  bathymetry  (especially  seabed  depth  under  the  large  Antarctic  ice  shelves)  and  rela¬ 
tively  sparse  and  poor  quality  data  for  model  validation  and  assimilation.  The  substantial  lack  of  data  in  polar 
regions  must  be  addressed.  We  suggest  that  more  bottom  pressure  gauges  should  be  deployed  at  specific 
locations  showing  enhanced  intermodel  differences.  As  an  example,  it  would  be  useful  to  deploy  bottom 
pressure  instruments  in  Baffin  Bay,  especially  in  Nares  Strait  connecting  Baffin  Bay  to  the  Arctic  Ocean.  New 
GPS  records  on  floating  Antarctic  ice  shelves  in  the  Weddell  Sea  have  recently  been  obtained  [King  et  al., 
2011a,  201 1  b];  comparable  records  over  the  Ross  Ice  Shelf  would  be  useful  to  resolve  uncertainties  around 
old  gravimeter  sites. 

We  expect  that  new  satellite  data  will  be  invaluable  additions  for  polar  regions.  ICESat  laser  altimeter  data, 
although  limited  by  its  short  time  series,  were  useful  in  the  southernmost  Antarctic  ice  shelves  [Ray,  2008; 
Padman  eta!.,  2008].  Existing  and  future  CryoSat-2,  FIY-2A,  Sentinel-3,  and  ICESat-2  data  will  lengthen  the 
time  series  and  should  considerably  strengthen  tidal  inversions.  Planned  wide-swath  altimeter  missions 
are  expected  to  overfly  some  polar  seas,  which  will  be  beneficial;  in  addition,  the  high  spatial  resolution  of 
such  altimetry  will  help  resolve  short-scale  tides  in  all  shallow  seas.  Notwithstanding  these  new  data,  satel¬ 
lites  over  polar  regions  are  unlikely  to  be  the  panacea  they  have  been  in  lower  latitudes,  for  several  reasons. 
(1)  All  high-inclination  satellites  have  inferior  tidal  aliasing  characteristics,  with  generally  one  or  more  con¬ 
stituents  aliased  to  very  long  periods.  (2)  Sea  ice  can  considerably  complicate  retrievals  of  sea  levels,  even  for 
smaller-footprint  altimeters.  (3)  Altimetry  over  floating  ice  shelves  will  remain  less  accurate  than  open-ocean 
altimetry  owing  to  unavoidable  changes  in  the  ice  and  snow  cover. 

Since  GRACE  data  appear  to  be  one  of  the  best  tools  for  identifying  high-latitude  tide  model  errors  (e.g., 
Figure  9),  these  data  in  principle  supply  valuable  information  for  constraining  models.  Some  early  efforts  at 
assimilating  such  data  have  begun  [Egbert  et  al.,  2009],  but  the  inherently  coarse  spatial  resolution  of  satel¬ 
lite  gravimetry  requires  that  any  assimilation  must  continue  to  give  heavy  weight  to  dynamical  constraints 
to  prevent  corrupting  shorter  scales.  Thus,  the  proper  exploitation  of  GRACE  data  will  necessitate  continued 
reliance  on  hydrodynamic  modeling,  with  its  associated  requirements  for  improved  dynamics  and  improved 
bathymetric  information. 

The  need  for  improved  bathymetric  data  sets  is  a  theme  that  runs  through  all  types  of  tide  modeling, 
whether  in  the  open  ocean,  shelf  areas,  or  high  latitudes.  It  has  long  been  known  that  bathymetry  is  espe¬ 
cially  important  for  shelf  tide  modeling,  but  the  required  accurate  data  sets  are  often  lacking.  The  estimation 
of  bathymetry  from  altimetry  data  works  much  less  well  in  shelf  regions  than  in  the  open  ocean  [Smith  and 
Sandwell,  1 997].  Thus,  our  knowledge  of  shelf  bathymetry  relies  on  acoustic  soundings,  which  for  many 
regions  of  the  world  are  not  in  the  public  domain.  In  polar  regions  updated  bathymetric  bathymetry  prod¬ 
ucts  are  being  developed,  using  a  combination  of  traditional  depth  sounder  and  swath  mapping  along  ship 
tracks  combined  with  interpretation  of  gravity  records  from  satellites  and  airborne  surveys  [Jakobsson  et 
al.,  2012;  Arndt  et  al.,  2013],  and  should  be  incorporated  into  global  tide  model  grids.  For  specific  regions 
of  large  tidal  error,  e.g.,  Larsen  C  Ice  Shelf  in  the  Weddell  Sea,  detailed  surveys  with  airborne  gravimetry 
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[Cochran  and  Bell,  201 2]  and  seismic  surveys  [Brisbourne  eta!.,  2014]  can  deliver  vastly  improved  water  col¬ 
umn  thickness  information.  Coordinated  international  efforts  to  improve  bathymetric  data  sets  remain 
crucial  for  tide  modeling. 

For  continued  advances  in  shelf  and  polar  regions  we  expect  data  assimilation  into  hydrodynamic  models 
to  play  a  larger  and  larger  role.  In  light  of  that,  it  is  encouraging  to  note  the  progress  that  has  been  achieved 
in  the  field  of  purely  hydrodynamic,  or  forward,  tide  modeling.  Such  models  have  advanced  in  recent  years 
as  a  result  of  more  powerful  computational  capabilities,  more  accurate  and  higher-resolution  bathymetric 
data,  and  better  physics  incorporating  findings  from  analyses  of  altimeter  data.  One  important  lesson  has 
been  the  need  for  a  gravity  wave  drag  mechanism  in  forward  models.  Improved  unconstrained  models  thus 
reflect  progress  in  understanding  tidal  dynamics.  For  the  M2  constituent,  totally  unconstrained  global  mod¬ 
els  are  now  about  as  accurate  as  the  circa  1 980  Schwiderski  solution  (see  Figure  1 2).  Notwithstanding  these 
advancements,  in  terms  of  accuracy  of  tidal  predictions,  unconstrained  models  are  not  competitive  with 
models  that  assimilate  altimetry. 

In  addition  to  improved  bathymetry,  tidal  hydrodynamic  modeling  is  intertwined  with  questions  concerning 
tidal  energy  dissipation  and  potential  interactions  with  the  oceanic  general  circulation.  The  various  param- 
eterizations  of  dissipation,  the  energy  loss  through  interactions  with  sea  ice  and  ice  shelves,  the  exchange 
of  energy  between  barotropic  and  baroclinic  motion,  and  the  physics  of  internal  wave  drag  remain  out¬ 
standing  topics  for  investigation.  Further,  we  need  improved  understanding  and  modeling  strategies  of  the 
seasonal  variations  of  barotropic  tides  induced  by  interactions  between  barotropic  tides  and  stratification  as 
well  as  seasonal  ice  coverage  [e.g.,  Rayetal.,  201 1;  Muller  et  al.,  2014;  Fok  etal.,  201 3]. 

Dynamically  consistent  data  assimilation  seems  a  key  way  forward  [e.g.,  Egbert,  1997;  Taguchi  et  al.,  2014], 
potentially  allowing  for  simultaneous  estimation  of  uncertain  model  parameters  (e.g.,  bottom  topography 
or  dissipation)  [cf.,  Losch  and  Wunsch,  2003].  Eventually,  the  validity  of  the  current  practice  of  separating 
tidal  motions  from  the  oceanic  general  circulation  will  start  to  break  down,  ultimately  requiring  ocean  and 
climate  models  to  have  tides  and  their  effects  directly  included.  Two  of  the  nonassimilative  forward  models 
we  tested  here  are  beginning  steps  in  that  direction. 

The  study  of  the  tides  began  long  ago  as  an  essential  part  of  nautical  life  [Cartwright,  1999].  The  subject 
led  to  breakthroughs  in  basic  geophysical  science  and  also  in  accurate  tidal  predictions  of  extraordinary 
practical  and  scientific  utility.  We  expect  that  future  advances  will  continue  to  yield  new  insights  into  the 
nature  of  the  ocean  while  continuing  the  long  heritage  of  remarkably  diverse  applications. 
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